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1. Introduction 

We have been engaged in developing a simple method to prepare by 
means of SiCl, a kind of silica-impregnated paper suitable for paper 
chromatography of phosphatides. Many comparative experiments were 
needed to ascertain the best method of preparing the paper, the compo- 
sition of a suitable mobile phase, and other characteristics of chromato- 
graphy with SiCl-treated paper!). For the purpose of comparative 
chromatography the two kinds of apparatus described in this paper were 
developed. Though especially adapted to the chromatography of phospha- 
tides on silica-impregnated papers with mixtures of di-isobutylketone, 
acetic acid and water (MARINETTI et al. [1]) as mobile phase, there are 
possibly other fields in paper chromatography in which the apparatus or 
some of its constructional features might also be useful. 
2. Requirements for a device, suitable for chromatography of phosphatides 


Previous experiments [2] on the chromatography of phosphatides on 
silica-impregnated paper, prepared according to Marinetti and colla- 
borators [1] with di-isobutylketone — acetic acid — H2O mixtures had shown 
that: 

a) The Re values strongly depend on the height to which the front of 
the mobile phase had ascended. 

b) The separation of the components of a mixture of phosphatides — the 
height of the front being equal—is improved with increasing volume 
of the vessel in which the chromatography is performed. 

c) Previous conditioning of the paper does not improve the results but 
on the contrary can have a deleterious effect. 


Our investigations on SiCl,-treated paper started with a conventional 
type of apparatus, shown in fig. 1. The mobile phase is present as a layer 
on the bottom of a cylindrical glass vessel of 570 ml capacity, closed by 
a ground stopper. To the lid a block of akulon (a kind of nylon) of the 
shape shown in the figure is fastened. The protruding ends of a stout 
—____—- mA 
1) To be published shortly in these Proceedings. e 
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platinum thread through the lower part of the block are bent into hooks. 
With this device two strips may be run simultaneously. Each strip has 
three perforations, the upper perforation, halfway the breadth of the 
paper, serves to hang the strip on the platinum hook. The strips remain 
stretched, and are kept separate, by means of a glass rod provided with 
four platinum hooks, which is fitted into the two lower perforations of 
each strip. 

With this apparatus we found that the statement (a) is valid also for 
the SiCl-treated paper. Rr-values, calculated in the usual way (distance 
displaced spot to starting point divided by distance front to starting 
point), are not constant when the travelling distance of the front is varied. 
Experiments indicated that another quotient is characteristic and con- 
stant, viz. the distance of a displaced spot to the immersion line of the 
strip in the mobile phase, divided by the distance of the front to this 
immersion line. This statement is valid, provided the front has ascended 


10 cm 


Fig. 1. Chromatographic apparatus of the conventianol type, allowing two strips 


to be run at the same time. 
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a sufficient distance above the Starting point. Remarkably enough, this 


was also found to apply to the chromatography of a single phosphatide 
with the same mobile phase on untreated paper. 


O 
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de 
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Fig. 2. Chromatograms obtained on non-impregnated Schleicher and Schiill paper 
with an apparatus of the type of fig. 1, in which lecithin has been applied at starting 
points a, b, ce and d. 


In fig. 2 we give the results of the chromatography of lecithin (phospha- 
tidyl choline) on untreated Schleicher and Schill paper No. 2043b. An 
apparatus similar to that represented by fig. 1 was used, the difference 
being only that the glass cylinder was longer, so that longer strips could 
be used. Two strips were suspended as in fig. 1, after about 5 mm? had 
been applied by means of a self-filling capillary micropipette of a 0.1 % 
solution of lecithin at the points marked a, b, c, and d1). The mobile 
phase was di-isobutylketone — acetic acid — water (40 : 25 : 5, by volume). 
The horizontal broken line drawn at the lower end of the two strips 
indicates the immersion line. After drying the strips were stained with 
the Acid Fuchsine—Uranyl—pH2 staining solution described else- 
where [3, 2]. After drying the chromatogram the red-coloured lecithin 
spots were outlined with a pencil, as well as the front which is made 
visible by U.V. light. 
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The lecithin spots which have migrated from the starting points a 
and b lie on one horizontal level. The spot migrated from the starting 
point ¢ has remained practically stationary on the same level (or moved 
to a very slightly higher level). Only the spot which migrated from starting 
point d lies distinctly higher. Thus in the cases a, b and practically ec, 
the position of the migrated spot is independent of the position of the 
starting point. In case d, where the front has not passed high enough 
above the starting point, the position of the migrated spot is no longer 
wholly independent of the position of the starting point. The values of 
Ry calculated in the usual way (distance “‘center’” of displaced spot to 
starting point divided by distance front to starting point) is not constant, 
clearly, but decreases greatly from a to d (a=0.70; b=0.63; e=0.51 and 
d= 0.34). On the contrary, the quotient: displaced spot to immersion line, 
divided by distance front to immersion line, is the same for a, b, and 
(practically) c (a=0.72; b=0.72; c=0.74 and d=0.77). 

In the communication directly following in this fascicule of the 
Proceedings, it will be shown that (provided the starting points lie 
sufficiently far below the front) the quotient: distance displaced spot to 
immersion line, divided by distance front to immersion line, is independent 
of the absolute value of the distance of the front to the immersion line. 

These results hold not only for untreated paper but also for SiCl- 
treated paper!). In constructing the chromatographic device we should, 
therefore, take care to fix the precise location of the immersion line of 
the strip. 

Point (b), mentioned at the beginning of this section, is illustrated by 
the following facts. 

We found that when SiCl,-treated paper strips of the same dimension 
are used in a much larger container the separation of the spots of the 
components is larger with the same ascension height of the front above 
the immersion line, just as in the case of silica-impregnated paper 2). 
Obviously (preferential) evaporation (probably of acetic acid) from 
the wetted strip into the available space plays a role. 

In constructing the chromatographic device we can take advantage of 
this evaporation effect. From the point of view of reproducibility this 
can best be realised by starting the chromatography in a vessel in which 
the air space does not contain vapours of the mobile phase. This can be 
achieved by bringing the mobile phase into the apparatus only just 
before starting chromatography. 

Besides, the mobile phase should have an as small as possible common 


) To be published later in these Proceedings. 

) The space in the container cannot be enlarged indefinitely as separation of 
the mobile phase in the strip will occur. This must be due to too great a loss of acetic 
acid to the surrounding space. The composition of the mobile phase in the strip 


changes as a consequence of loss of acetic acid in such a way that ultimately sepa- 
ration into two layers will occur. 
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surface with the air. This leads to the characteristic feature of the types 
of apparatus described below, namely a horizontal tube (inner diameter 
about 7 mm) provided with a slit (about 1-1.5 mm wide) which serves as 
container for the mobile phase. 

The paper strip or sheet dips with its lower end through this slit into 
the interior of the slit-tube. This solution of the feeding problem has 
two additional advantages: 

1. We may consider the entrance of the slit to indicate the immersion 
line on the paper when, after filling the slit-tube with the mobile phase, 
chromatography begins. The position of the immersion line should be 
known (see above). We draw therefore with a pencil a horizontal line on 
the paper strip or sheet, a small distance from the lower margin (5 mm), 
which represents the future immersion line. Then the suspended strip or 
sheet is let down into the slit until the pencil line coincides with the slit 
entrance. For reasons to be discussed later, the lower end of the strip 
should not reach to the bottom of the slit-tube, only to its axis. 

2. The small bore (7 mm) of the slit-tube and the bent filling tube 
(sealed to it midway at the opposite side of the slit) means that chromato- 
grams are run with a relatively small quantity of mobile phase. This 
allowed an economic use of mobile phase. We adopted the rule to run 
each chromatogram with a fresh quantity of our stock of mobile phase. 


3. Slit-feeding apparatus for relatively small-sized strips 


General 

For much comparative research work the apparatus given in fig. 38, 
A and B has been used. The front is allowed to ascend only 17.5 above 
the immersion line. This saves much time because the ascension time 
increases quadratically with the height of ascension. With SCHLEICHER 
and ScHiiLL paper no. 20436, untreated or treated with SiCly, the 
ascension time at 20° remains below 3 hours. 

The cylindrical vessel A can be closed with a loosely fitting paraffined 
cork B which has three holes, one for the filling tube C carrying the slit 
tube D, the second for an easily sliding glass rod E ending in a hook F 
and third for a short glass tube G, which serves as an outlet for air when 
the mobile phase is pipetted in. After filling it is covered, e.g. with a 
small glass plate. The capacity of the glass cylinder is 660 ml. 


Hanging strips in the right position 

We use strips 20.5 cm long and 4.5 em wide (fig. 4A). The strips have 
three pencil-lines, the future immersion line i, 5 mm above the lower 
margin, a starting line, 25 mm higher and a line, f, 17.5 cm above the 
immersion line. To line f the front will be allowed to ascend. It is con- 
venient to draw a fourth line on the strip 12.5 cm above the immersion 
line which serves to calculate approximately the moment when the front 
will reach the drawn frontline and the strip must be taken out. The time 


or) 


Fig. 3. Slit-feeding apparatus for relatively small-sized strips (4.5 x 20.5 em) 
(see text). 


which elapses for the ascension of the front to the line drawn 12.5 em 
above the immersion line is noted. It will then take practically an equal 
time for the front to reach the 17.5 cm line }). 

After having applied the spots on the starting points (crosses on line 8), 
the cork is taken out of the glass cylinder. The strip is hung with its 
perforation on glasshook F and the lower end is put through the slit 


1) The height of ascension in our relatively short strip is practically proportional to 
the square root of the ascension time. When a time ¢ is needed for an ascension height 
of 12.5 em, it will take very nearly a time 2t (reckoned from the beginning) for 
reaching 17.5 cm (because 17.5/12.5 = 1.40, which approximates to /2 = 1.41). 


Thus the time from start to 12.5 line is practically equal to the time from 12.5 em 
line to 17.5 em line, 
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(47 mm long) of the slit-tube D. The immersion line i drawn on the paper 
should coincide with the entrance of the slit tube. When the strip does 
not have its correct position it may be lowered or drawn up by means 
of the easily sliding rod E. 

The cork is replaced on the glass cylinder and eventual corrections in 
the position of the strip are made by means of rod E (up- or downwards 
for the right position of the immersion line. turning round its axis for 
correcting an eventually not parallel position of the strip to the walls of 
the glass cylinder). 


Fig. 4A. Strip to be used in the slit-feeding apparatus of fig. 3. 
Fig. 4B. Chromatogram obtained on SiCl4-treated paper with the apparatus of 
fig. 3. It is discussed at the end of section 4. 


General precautions and provisions 

The chromatography is performed in a dark room, the temperature of 
which is kept constant at 20°. If radiation from one side reaches the 
apparatus, the ascending frontline will not be horizontal; this should be 
avoided. Behind the chromatographic apparatus a small incandescent 
lamp (4 V.) may be switched on, to control the manner in which at the 
start of the chromatography the front passes the slit (see below) and to 
check from time to time the position of the front. 


Filling the slit-tube 
We now proceed to fill the slit-tube. For our apparatus some 3 ml of 
the mobile phase are needed. But this must not be pipetted at once into 
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the funnel of the filling tube. We have already mentioned in the preceding 
section sub 1) that the lower margin of the paper strip must not rest on 
the bottom of the feeding tube but should reach about the axis of this 
tube. no farther. Filling proceeds by slowly adding with a pipette a 
first portion of the mobile phase, which thus may have time to spread 
over the whole bottom of the slit-tube without as yet wetting the paper. 
Then, less slowly, a further portion is added, still insufficient to fill the 
slit-tube, but enough to wet the whole lower rim of the paperstrip. We 
then wait until the front has passed outside. In this way it is easy to 
achieve that the front, when it becomes visible above the slit, is reasonably 
horizontal. When the second portion added is too large, mobile phase 
may locally enter the slit; an irregular front will be the result. Though 
often the front becomes nearly horizontal again before reaching the 
spots on the starting line, it is better to avoid an irregular front at the 
very beginning. 

When the frontline has thus passed horizontally out of the slit, the 
slit-tube is filled further. One should avoid adding too much mobile 
phase, so that it streams out of the slit. In particular one should not 
try to get rid of an air bubble remaining in the slit-tube by adding more 
mobile phase. By capillary force the slit will then be sealed already with 
the mobile phase, so the paper is wetted everywhere in the slit. When 
more mobile phase is added the danger exists that mobile phase makes 
its exit through the slit. 


Cleaning after use 


After having finished chromatography, the apparatus must be cleaned. 
The cork is taken out, the remainder of the mobile phase is poured out 
and discarded. We then pour alcohol into the funnel, so that it flows 
out of the slit. When not needed directly for another run, the unit is 
hung upside down; the next morning it will be perfectly dry. When needed 
immediately for a new run, the unit is rinsed with acetone and dried with 
a Fohn. 

The cylindrical vessel is put down obliquely on its brim, in order that 
the heavy vapours may fall out. As a rule, the following morning it is 
free from vapours. When directly needed for another chromatographic 
run, the cylinder is freed from vapours with a Féhn. 


Use of more apparatus at the same time 


We have four apparatus of this type. Their dimensions should be as 
closely equal as possible. In particular the slit-tubes should be adjusted 
so that their distance to the bottom of the glass cylinder is equal. Despite 
such precautions minor differences may influence the degree of dispersion 
of the spots on the chromatogram. It is thus advisable to check this by 


running with them strips with the same phosphatide mixture brought up 
on the starting points, 
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4. Slit-feeding apparatus for paper sheets 
General 


This much larger apparatus, shown by fig. 5, A and B, in which paper 
sheets 17 cm wide and 27 em long are used, allows us to compare a larger 
number of spots in one run. The front can ascend more than 17.5 cm 
above the immersion line the latter drawn again 5 mm above the lower rim 
of the sheet. 

The volume of the glass container A (5,91) is much larger than that of 
the apparatus of fig. 3. This is why, with equal ascension height of the 
front, chromatograms are obtained in which the spots corresponding to 
components of a phosphatide mixture lie farther apart. The container 
has a broad rim, which is ground flat, so that no visible space remains 
when the glass plate B (of plate glass) serves as a cover. 

The U-shaped frame CDE. attached to the glass plate by means of 
corks which are inserted in appropriate holes, carries the slit-tube F. 
The left half CD of this frame consists of a glass tube entering into the 
slit-tube, which serves as filling tube. The right half DE is a solid glass rod 


10 cm 


Slit-feeding apparatus for paper-sheets 17 x 27 cm. Because of the much 


Fig. 5. a 
larger space of the container resolution is better than with the apparatus of fig. 3. 
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of the same outer diameter as the tube CD. The corks inserted into the 
holes of the glass plate allow the tubes C and E to be slid up- or down- 
ward, but not too easily. Care is taken that C and E are adjusted in the 
holes of the corks in such a way, that the slit-tube is exactly parallel to 
the upper surface of the glassplate. 

When the glassplate is placed as a lid on the container A, we can 
control whether the slit really occupies a horizontal position, by placing 
a water level on the lid parallel to the slit-tube. When necessary the 
container is lifted on one side by placing something underneath until 
the water level on the lid indicates a horizontal position of the slit. 
Furthermore the glass plate carries two glass rods G with hooks H by 
which the paper sheet is suspended. The glass rods can slide easily in the 
holes of the corks to bring the paper sheet in the position desired. 


Technique 


The technique in working with this larger apparatus is in principle the 
same as the one applied for the small apparatus described in the preceding 
section and need not be repeated here. The paper sheet, provided with 
immersion line (5 mm higher than the lower margin of the paper), starting 
line, future front line and possibly a fourth line indicating about half of 
the total ascension time 1s now suspended by two perforations on the 
hooks H. Manipulation with both rods G is now generally necessary to 
bring the paper sheet in the position desired. The immersion line should 
coincide with the entrance of the slit, which can be adjusted by moving G 
up- or downwards; the left and right margin of the paper should not 
touch the ends of the slit which can be adjusted by turning the axis. 
The sheet is 170 mm wide, and the slit has a length of about 172 mm. 

To manipulate the unhandy large glass plate with attached frame, the 
stand pictured in fig. 6 is used. It consists of a wooden ground plate A, 
and two broad (13 cm) aluminum strips, B, on the back reinforced by 
two crossed aluminum strips C., 

The broad strips are bent at their upper end in the shape as given in 
the figure. To bring the paper sheet in place the glass plate is laid down 
on the stand as shown. 

The plate is now placed on the container. We take care that the glass 
frame takes in a symmetrical position (by viewing in three directions). 
We now control with the water level the horizontal position of the slit 
and if necessary make a last correction of the position of the paper sheet. 
We then proceed to fill in the necessary amount of the mobile phase (in 
our apparatus we need in total about 7 ml) following the directions given 
in the preceding section. When everything is done properly the front 
rises as a horizontal line out of the slit. 

In cleaning the apparatus the glassplate with frame is held upside 
down, to remove most of the remaining mobile phase (which is discarded), 
Then the plate is placed on the stand, in the position as shown in the 
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figure, and the filling tube and slit-tube are rinsed thoroughly by pouring 
alcohol into C by means of a funnel. A tray placed on the wooden Sieee 
plate catches the alcohol streaming out of the slit. The glass plate with 
attached frame is now placed upside down on the frame to dry overnight. 

When we want to continue chromatography immediately, rinsing with 
alcohol is followed by rinsing with acetone and drying with a Fohn in 
upside down position on the stand. 


LL} 
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Fig. 6. Stand for conveniently manipulating the glass plate with attached frame 
Olio o. 


Check on functioning of the apparatus 

In fig. 7 we give two chromatograms (stained with Acid Fuchsin and 
Uranylnitrate at pH=2) on SiCl,—treated ScuiuercHeR and ScHULL 
2043 b. The mobile phase was di-isobutylketone—acetic acid — water 
(50 : 25 : 5 by volume). The temperature was 20°. At the starting points 
5 mm3 quantities of a 1 % or 2 % phosphatide solution were applied by 
means of a self-filling capillary pipette. The phosphatide solution contained 
a mixture of lecithin (phosphatidyl choline), cephalin (phosphatidy] 
ethanolamine), their lysoproducts, and sphingomyelin. The solution was 
obtained by mixing column fractions in suitable proportions, and it was 
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Fig. 7. Chromatograms on SiCl4-treated Schleicher and Schiill paper (20436), of 

a mixture of phosphatides (C = cephalin; L = lecithin; LC = lysocephalin; LL = 

lysolecithin and S = sphingomyelin). No separation is obtained of LC and 8, The 

numbers | and 2 refer to the total concentration (°,) of the mixed solution of which 
5 mm® was brought up on the starting points. 


generally used as a reference mixture on paper chromatograms !). On 
silica-impregnated paper four spots are obtained in the order: (front to 
starting line) cephalin (C), lecithin (L), lysocephalin plus sphingomyelin 
(LC+$) coinciding, and lysolecithin (LL). The sequence of spots on 
SiCl,-treated paper is the same. 

The chromatograms of fig. 7 differ in degree of resolution, which is 
clearly greater on the left chromatogram. Why resolution should differ 
can only be discussed in a future communication, dealing with preparation 
and properties of SiCly-treated papers. At this point the two chromato- 
grams are given to illustrate the functioning of the apparatus. If the 
apparatus works ideally, and if the paper has over its whole area the 
same properties, it is to be expected: 

a) that the front should be horizontal, and 


b) that spots belonging to the same component of thes mixture lie on 
one horizontal line. 


') We thank Dr. G. J. M. Hoocuwtnket for putting this reference mixture at 
our clisposal. 
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The chromatograms of fig. 6 show that these expectations are reasonably 
fulfilled. The distance between the starting points differs in the two 
cases; on the left chromatogram this distance is 2 em, on the right 1.5 em. 
Even with 1.5 em distance, where 11 spots were applied, no disturbance 
occurred. The conclusion is that the apparatus functions adequately. 
Parenthetically we may also conclude that SiCly-treatment has apparently 
resulted in a homogeneous paper over its full width. Attention should 
further be drawn to the fact that with the slit-feeding apparatus used 
here (fig. 5) the resolution of the components of a phosphatide mixture 
is better than in the case of the small slit-feeding apparatus (fig. 3) with 
much smaller air space. Here the chromatogram of fig. 6B should be 
compared with that of fig. 4B, since in both cases the same kind of SiCly- 
treated paper (S & S 2043 b) was used, and the front in both cases had 
ascended the same distance (17.5 cm) from the immersion line. 


Summary 


1. In constructing a device suited for reproducible and comparable 
chromatography of phosphatides (silica-impregnated paper; di-iso- 
butylketone — acetic acid —H20 as mobile phase) one must consider with 
the following facts: 

a) Re values strongly depend on the ascension height of the front. 

b) The place of the immersion line on the chromatogram is (in practical 
cases) of much more importance for characterizing chromatographic 
displacement than the place of the starting spot. 

c) Previous conditioning of the paper is not required. 

d) Separation of the spots of a mixture of phosphatides is better when 
working in a vessel with larger air space. 

2. The requirements for a chromatographic device are, accordingly: 
a) The immersion line on the strip should be known. 

b) The state of the apparatus at the beginning of each run should be 
exactly reproducible. 

These requirements are fulfilled in the type of apparatus described. 
Its characteristic feature is that the strip or sheet dips through a slit into 
a feeding tube, (whereby the immersion line is defined), which is provided 
at the beginning of the run with the mobile phase through a filling tube. 
After each run the remaining mobile phase is discarded, and the apparatus 
cleaned, so at the beginning of a new run the apparatus is in the same 
state as at the beginning of the previous run. 

3. Two types of slit-feeding apparatus are described: 

a) a type for small strips (4.5 x 20.5 cm), allowing an ascension height 
of the front of 17.5 cm above the immersion line. 

b) a type for paper sheets (17x 27 cm), allowing, if necessary, a higher 
ascension of the tront. 

4. Directions are given concerning the use of the apparatus, and 
precautions necessary to obtain reproducible chromatograms are indicated. 
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Some examples are given pertaining to the chromatography of 


phosphatides demonstrating : 


a) 


b) 


that the apparatus for paper sheets—despite its large dimensions — 
functions satisfactorily. The front is horizontal, and spots corresponding 


to one component lie on one line. 
that this apparatus, because of its much greater air space, gives a better 
resolution of the spots than the apparatus for small strips. 


Department of Medical Chemistry, 
University of Leyden 
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BIOCHEMISTRY 


CHROMATOGRAPHY OF LECITHIN ON NON-IMPREGNATED 
PAPER BY MEANS OF A SLIT-FEEDING APPARATUS 


BY 


H. G. BUNGENBERG DE JONG anp J. To. HOOGEVEEN 


(Communicated at the meeting of September 26, 1959) 


Introduction 


Certain types of slit-feeding apparatus have been developed (see 
preceding communication in this fascicule) to meet the requirements for 
reproducible chromatography of phosphatides on impregnated paper [1]. 

The present investigation provides an example of the suitability of 
the large type of slit-feeding apparatus for systematic investigations. 
As shown in the title non-impregnated paper is used here. It would seem 
that the results of the present investigation have no significance for 
chromatography of phosphatides on silica-impregnated paper. In later 
publications on the use of SiCly-treated paper in chromatography of 
phosphatides, the results will, however, prove useful and enable us to 
draw certain general conclusions. 


2. Hapervmental 

Five sheets, 17 x 28 cm, were cut in the same direction out of SCHLEICHER 
and ScHULL paper no. 2043, the long side of the sheet lying in the direction 
of the smaller ascension rate. With a pencil were drawn: the immersion 
lines (5 mm above the lower margin of the paper), the “front lines” (at 
a distance from the lower margin of the paper of 5; 10; 15; 20 and 25 cm 
respectively), and lines proceeding upward to the right with nine starting 
points on equal horizontal distances. 

On each starting point 5 mm? of a 0.1 % solution of lecithin was 
applied. The five chromatograms were run in succession, using di-isobutyl- 
ketone —acetic acid—H2O (50: 25:5, by volume) as mobile phase. 
They were stained with Acid Fuchsin—Uranyl nitrate at pH=2 [2, 3]. 
After renewed drying, the red spots were outlined by means of a pencil 
At the same time the position of the front lines was indicated. 


3. Discussion of the results 

The five chromatograms obtained are reproduced in the figures 1-5. 
The migrated lecithin spots do not lie on a straight line drawn through 
the intersection point of the frontline with the line ascending through 
the starting points. As any straight line through the above intersection 
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Fig. | Fig. 2 


Figs. 1 and 2. Chromatograms of lecithin on non-impregnated paper, provided 
with nine starting points on sloping lines. The front has risen to approximately, 
4.5, and 9.5 em above the immersion line. 


point (and proceeding between front and line through the starting points) 
represents a definite Rf value, it is evident that R/ values (characteristic 
in the case of ideal partition chromatography) are far from constant. 
The curves through the transported spots tend to reach a horizontal line 
on the left side of the chromatograms. It follows that on the left side the 
position of the transported spots tends to be independent of the position 
of the starting points. 

This may be an example of chromatography by means of a gradient 
in the composition of the mobile phase, which develops in the paper 
during the upward movement of the mobile phase. The horizontal lines 
which the spots tend to approximate asymptotically on each of the 
chromatograms (figs. 1-5) may indicate a certain characteristic compo- 
sition within the gradient of the mobile phase. For this composition the 
upward velocity of the lecithin spot, divided by the velocity of the front, 
probably reaches its largest value. 

Assuming that the two spots farthest to the left on each chromatogram 
lie already on this asymptote, we may investigate the position of this 
asymptote as a function of the position of the front. Table I gives in 
column 2 under the heading “asymptote” the mean distance of the 
centres of the two left spots to the immersion line. Similarly column 3 
under the heading “front” gives the mean distance of the front (above 
the two spots) to the immersion line. These distances may be designed 
by the symbols a; and f;. The suffix i indicates that the distance of the 
asymptote (a) and of the front (/) is measured here from the immersion 
line (7). In fig. 6 the values a; (column 2 of the Table) have been plotted 
against the values /; (column 3). The straight line drawn through the 
origin reasonably fits in with the experimental points, thus indicating 
that a is proportional to fj. 

The values of the quotient a;/f; in column 3 of Table I differ very little. 
The differences lie within the experimental error. Considering the quotient 
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a;/f; as constant we may introduce for it a symbol Rfi, which is defined: 


Rina a, _ distance of asymptote to immersion line 
me hi distance of front to immersion line 


For the apparatus used here, the front height evidently is again pro- 
portional to the square root of ascension time (columns 4 and 5, Table I). 


TABLE I 
Distance from immersion line to horizontal line through the two left spots 
(“asymptote”) and to the front (‘‘front’’) 


a 


Chromatogram | ““Asymptote” | “Front” : Ascension time | I 
= E : | Rfi = ai! fi - . 3 fil Vt 
figure im mm , In mm In minutes 
1 29 | 46 0.63 | 12 13.3 
2 62 96 0.65 47 14 
3 95 145 0.66 110 | 13.8 
4 127 194 0.65 193 | 14 
5 165 245 | 0.67 | 332 | 18.4 
| mean 0.65 | | mean 13.7 


Above we made the assumption that on each chromatogram the position 
of the asymptote (a;) denotes that particular composition of the mobile 
phase at which lecithin migrates upward at maximum speed, compared 
with the velocity of the front. The constancy of Rf; indicates that the 
above particular composition within the gradient is always found at 
the same fraction of the height of ascension of the front (f;), and we 
suggest that for other compositions within the gradient of the mobile 
phase the same may hold. We will return to this question in a later 
communication dealing with chromatography of phosphatides on SiCl,- 
treated paper. 
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Fig. 6. ‘“‘Asymptote” (aj) as a function of “front” (fi). 


We will finish this communication with two remarks: 

1. The motive for the present systematic investigation is derived 
from a single experiment on the chromatography of lecithin on non- 
impregnated paper in the preceding communication, fig. 3, using as 
chromatographic apparatus a closed tall glass cylinder with a layer of 
mobile phase at the bottom. It was found that with sufficient height of 
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the front the position of the migrated lecithin spot is independent of the 
position of the starting point. Here the same conclusion is arrived at, 
witness figs. 1-5. We may conclude that the peculiarities which have been 
found in the present investigation were not introduced by the slit-feeding 
apparatus itself. 

2. It may be expected that the value of Rf;, found here to be 0.65 for 
lecithin, will depend on the dimensions of the slit-feeding apparatus used. 
In particular the volume of the air space of the container will influence 
the position and steepness of the gradient, since evaporation from the 
wetted paper will be different. This may thus lead to another value of 
ff;. So for systematic investigations in which chromatograms must be 
compared, it is safest to use one and the same apparatus. 


Summary 


1. The chromatography of lecithin on non-impregnated paper with 
di-isobutylketone — acetic acid—HzO (50: 25:5 by volume) has been 
studied by means of a large slit-feeding apparatus. Five chromatograms 
were run, in which the front was allowed to ascend to different distances 
above the immersion line. On each chromatogram the lecithin was applied 
on nine starting points, lying on equal distances on a sloping line inter- 
secting the front on the paper sheet. 

2. The migrated lecithin spots on each chromatogram are situated on 
a curve, in such a way that the lowest starting points reach a practically 
horizontal line (asymptote). 

3. Since a straight line cannot be drawn through the intersection point 
of front and sloping line and through the centres of all spots, Rf values 
(in the ordinary sense) are not constant. The Rf value therefore is not a 
suitable number for characterising a particular substance. 

4. The quotient of distance a; of the asymptote mentioned sub 2. to 
the immersion line, divided by the distance f; of the front to the immersion 
line, is independent of the ascension height of the front. It is proposed 
to denote this quotient with the symbol Ff;,, thus Rf;=ai/fi. 

5. The results suggest that a gradient is present in the mobile phase. 
The asymptote may indicate a particular composition of the mobile phase 
at which the upward migration of lecithin reaches its maximum velocity, 
relative to the velocity of the front. 

6. The experimental fact that Rf; is constant may mean that this 
particular composition is always situated at a constant fraction of the 
distance from front to immersion line. 
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PALEONTOLOGY 


REMARKS ON A FOSSIL HUMAN MOLAR FROM OLDUVAT, 
EAST AFRICA 


BY 


G. H. R. VON KOENIGSWALD 


(Communicated at the meeting of October 31, 1959) 


One of the most important prehistoric sites in Africa is Olduvai Gorge 
in Tanganyika. Here we find, beautifully exposed, extensive series of old 
lake and river deposits, containing not only a rich fauna, but also abundant 
industry. In (at least) ten levels we can follow the transition from the 
Oldowan pebble culture through primitive and advanced Chellean into 
the Acheulean. With great care and success Dr. L. 8. B. Leakey, assisted 
by his wife, M. Leakey, has directed excavations and, collecting at 
Olduvai, brought together an enormous wealth of material. 

While stone implements are plenty, human skeletal material is extremely 
rare. In 1955 two human teeth were discovered in a Chellean I living site 
at BK IT in the lowest level of bed II at Olduvai. As no additional material 
came to light, they have only recently been described by LEakKgEy (Nature, 
1958). One is a small canine, the other a large molar of complicated 
pattern. 

The interpretation of isolated human teeth is not always easy. While 
Leakey “after very lengthy study and comparison” is satisfied that the 
molar is a left lower second deciduous molar’ of Homo, he freely admits 
that “there are also those who maintain that we are dealing with an upper 
permanent molar of australopithecine type”’. 

During the Pan African Congress in Livingstone the author had the 
privilege of seeing the original material and his first impression has also 
been that it must be an upper molar. Engaged in a study on the dentition 
of Karly Man, he was puzzled by the photographs published and struck 
by the remark that such an early human molar should have a cusp pattern 
“diametrically opposite to the so-called’? dryopithecus pattern.’’ We have 
to thank the British Museum (N.H.), which upon request kindly put an 
excellent cast at his disposal, which forms the basis for this short study. 


The roots. 


Fortunately, an essential part of the roots has been preserved in the 
fossil molar. As the roots of upper and lower molars exhibit principal 
differences as well in number as in position, we might study this part of 
the tooth first. 


To quote from Zrisz and NuCKHOLLS (1949) the root of the first permanent 
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lower molar of man—the deciduous teeth principally show the same 
arrangement —“‘presents a trunk which bifurcates very close to the 
cervical line to form two branches which are named for the positions they 
occupy, the mesial and distal terminal roots or branches’. The root of 
the first upper molar ‘‘consists of a trunk which trifureates a short distance 
from the cervical line to form three branches, or terminal roots, which 
are named for the positions they occupy, as the lingual, the mesio-buccal 
and the disto-buccal branches’’. To illustrate these conditions, we have 
here given the outlines of a last upper (figs. 1, a-c) and a lower (figs 2, a—-c) 
deciduous molar after MARSEILLIER (1937). 

A glance at the molar from Olduvai (figs. A, B) reveals that there are 
three roots and that in number, arrangement and direction they correspond 
exactly to the roots of an wpper molar. The two smaller and nearly straight 
roots mark the buccal side (fig. A); the large root, set in an angle, must 
indicate the lingual side (fig. B). Such divergent roots, which, when the 
crown is viewed from above, stick out lingually, occur in anthropoids and 
in Karly Man. WEIDENREICH has figured two comparable molars of Peking 
Man (1937, figs. 109 and 111), in which the lingual root shows nearly the 
same inclination; one of the molars is copied here in fig. 3. From a drawing 
by SELENK4 for orang (1898, fig. 105) it is evident that this angle is greatest 
in the first, less in the second and absent in the third molar; both molars 
of Peking Man are first ones. This might be taken as an indication that 
our molar in question must be regarded as a first molar. As the lingual 
root is slightly directed towards the rear (compare our fig. la), this root 
furthermore indicates a left molar). 

Remains the last question: is this a deciduous or a permanent molar? 
With one isolated tooth, this is difficult to answer. Are the roots really so 
delicate? As differences between the roots we found: that in permanent 
molars the roots are longer, and that they are less arched; in deciduous 
molars they are more spaced to leave room for the premolar. We have 
sketched the situation in fig. 4. From this it is most probable that our 
molar is a permanent one; the last reason for such an interpretation is 
that we are dealing with a very large tooth; if it were a deciduous molar, 
a permanent would be still larger, which seems improbable in a hominid, 
even with regard to the robust dentitions of the Australopithecinae. 

A study of the roots leads us to the conclusion that the molar from 
Olduvai is a left wpper first permanent molar. 


The crown. 

This is confirmed by the crown itself. The occlusal surface shows a 
normal pattern, except for a “supernumerary cusp” (LEAKEY), — fig. C2, a 
—which blocks the outside of the valley between protocone (pr.) and 
hypocone (hy.). Paracone (pa.) and metacone (mt.) are of about equal 
size; the latter is connected with the protocone by a slightly interrupted 
crista obliqua (co.), as is shown in fig. C. 


Bigs, AY B; 6. 


Figs. 1, 
ioe. 
Fie. 4. 


All tigures in 
Gg 4. i a% 2/1 nat. size 


First upper molar of Australopithecus from Olduvai. A buceal, ; 
B, mesial view (to demonstrate the roots); C, occlusal view (1) 


pretation (2) of the pattern. 
2. Last left deciduous molars of man (to demonstrate the position of 


the roots). 1) upper, 2) lower molar, a lingual, 6 buceal, ¢ mesial view. 
(Redrawn after MARSEILLTER 1937, reversed). 
Sinanthropus pekinensis, first upper molar (nr. 95) 
WEIDENREICH 1937), 


with inter- 


in mesial view. (After 
Differences in size and roots of the first permanent (full line) and last 


deciduous molar of man. Buccal view. Drawing after data from ZEISZ 


and Nucknonyis 1949, 
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The anterior fovea (fa.) is wide. The anterior border is slightly beaded, 
with a more marked cusp at the lingual border, which sometimes is called 
a “protoconulus” (pr.). Between this cusp and the protocone there is a 
slight interruption, which at the lower border ends in a “Carabelli pit” 
(cp.), which is very typical for the Hominidae. The posterior fovea (f.p.) 
is elongated too, there is a fissure in the posterior border, with secondary 
cusps (s.) on both sides. 

There is nothing unusual in the Olduvai molar, except its relative 
dimensions. In the upper molars of the Hominidae the buccolingual 
breadth generally exceeds the mesiodistal length; here we see the reverse. 
This is caused by the well-developed fovea, and it is probably this 
lengthening of the tooth in mesiodistal direction, which has given the 
wrong impression that the tooth might be a lower molar. 

As has already been mentioned, in upper molars the breadth exceeds 
the length. But there are exceptions, and it is for our case important to 
see that we do not find these in Peking Man or Pithecanthropus, but in the 
Australopithecinae of South Africa, which in size and crown pattern show 
the closest resemblance to the Olduvai tooth. 

If we compare our molar with the first molars of Australopithecus 
crassidens from Swartkrans, its buccolingual breadth falls well into the 
limits of this tooth, but the mesiodistal length is greater. The difference 
is only 0.5 mm and not significant. Among the first molars from Swartkrans 
there are also four (Sk. 27, 47, 55, 839) in which the length exceeds the 
breadth up to a maximum of 0.8 mm; in our case it might be 1.0 mm. 


Length and breadth of upper molars 


d4 M2 
il. br. 1. br. 
Australopithecus africanus 10.1 11.0 12.8 14.0 
Australopithecus crassidens (av.) 10.6 IEG 13.8 14.5 
Sk. 838 11.0 12.1 
Sk. 89 14.6 15.0 
Sk. 27 14.0 Tee) 
Olduvai 15.0 14.0 
Pithecanthropus modjokertensis 12.0 14.0 
Sinanthropus pekinensis (av.) 10.9 12.5 
nr. 33 sit 13.4 
Homo neanderthalensis (Krapina) 8.5-10.0 10.0-11.38 10.0-11.8 11.5—-12.0 
8.2 10.0 10.5 11.8 mm. 


Homo sapiens (av.) 


Our tooth seems not to have been used yet. All the cusps are clear; 
the surface is not completely smooth, but there are not many wrinkles. 

Besides the molar, and apparently belonging to the same individual, 
a small canine has been found, which indeed might belong to the deciduous 
dentition. The tooth is worn, the rootsection is elliptic, and there is nothing 
against regarding the canine as an upper tooth. In size it is comparable 
to the same tooth of Australopithecus africanus; both dentitions might 
have been in the same stage of eruption and of the same individual age. 


24 


Size of deciduous upper canines 


length breadth 
Australopithecus africanus 6.8 5.8-6.0 
Olduvai 7.0 6.1 mm. 


So our deciduous canine has the same size as in modern Man; as, how- 
ever, the molars are much larger, it is relatively smaller. 

We have tricd to estimate the relative size of the canine in comparing 
its mesio-distal length with the same length of the last deciduous molar 
and the first permanent molar respectively. For Olduvai we get an index 
of 214; for Australopithecus africanus 148 and 187, for Homo sapiens 117 
and 153. The indices are higher for the permanent molar than for the 
deciduous one, and they are also higher for Australopithecus than for 
Homo. Consequently we see a confirmation of our opinion that we are 
dealing with a first permanent molar rather than with a deciduous tooth, 
and with an Australopithecine and not with Homo. 


Conclusions 

Are we dealing with the teeth of Homo or of Australopithecus? Or to 
put the question more concrete (as in our opinion both belong to separate 
lines of evolution), to an early Homo or a late Australopithecus? 

Leakey already tried to compare the material under description with 
Atlanthropus from Ternifine, the maker of early handaxes in North Africa. 
Unfortunately no upper molars are known, but there can be little doubt, 
as can be concluded from the lower jaws, that the dentition must be about 
the same as in Sinanthropus. Now if we make a closer comparison between 
the molar from Olduvai and those of Sinanthropus (WEIDENREICH 1937) 
and the Australopithecinae (RoBINsON 1956), size and morphological 
affinities clearly favour the Australopithecinae. In pattern our molar is 
nearly identical with a second molar of Australopithecus crassidens from 
Swartkrans, Sk. 831 (Roprnson, fig. 25d), which shows a Carabelli pit 
and the same accessory cusps in the posterior fovea. Also the size relations 
of the deciduous canine point in the direction of the Australopithecine. 

We therefore come to the conclusion that the upper molar from Olduvai 
must belong to an Australopithecus, probably related to A. crassidens. 
The presence of an Australopithecus together with Chellean Man is certainly 
a great surprise, but it is in accordance with the palaeontological and 
geological results of Ewmr, OakLEy, Howexu and others. The classical 
south african Australopithecinae are of Lower Pleistocene age; the 
Australopithecine from Olduvai is most probably the latest survivor of 
that group in Central Africa. 


Post scriptum: 


After this manuscript had been prepared, Dr. LzEAKEY announced the 
discovery of a complete skull from the upper part of bed T in Olduvai. 
(Illustrated London News, September 12 and September 19, 1959; Nature, 
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August, 15, 1959). Zinjanthropus boisei LuaKkny, according to its dentition 
with enormous molars and very reduced canines, in the possession of a 
sagital crest, it is the most specialised member of the Australopithecinae, 
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GEOPHYSICS 


THE DIFFERENCE OF THE TECTONIC DEVELOPMENT ON 
THE EAST AND THE WEST SIDE OF THE PACIFIC 


BY 


F. A. VENING MEINESZ 


(Communicated at the meeting of October 31, 1959) 


Before undertaking an attempt to explain the great difference of the 
tectonic development on the east and the west side of the Pacific, where 
we find the Rocky Mountains-Andes orogenesis on one side, and the 
island-ares to the east of Asia and Australia on the other, we must adopt 
a hypothesis about the cause of the great tectonic deformations of the 
earth’s crust in general. We shall assume as such the hypothesis of systems 
of convection-currents in the mantle of the earth, which exert strong drag 
forces on the crust. We shall not here make mention of the great many 
arguments in favour of this hypothesis, for this purpose we may refer to 
previous papers by the writer!); they make it highly probable, if not 
certain. 

Considering this hypothesis we must realize that it is likely that the 
mantle over its whole depth of about 2900 km is crystalline. Currents 
in the mantle must, therefore, have the character of deformations in solid 
crystalline matter, and so these deformations must be of a plastic kind. 
This involves an elastic limit which the stress-deviator has to exceed 
before permanent deformation and flow is possible, and this has two 
consequences, 

In the first place a rising current can in the surface layer of the mantle 
not flow out to all sides, as usually a hydrodynamic current in a viscous 
Newton fluid must do, but it can only do so to one side. This is caused 
by the fact that once the elastic limit is overcome in one direction, flow 
sets in towards this direction and this renders it unlikely that in another 
direction the elastic limit is attained. This onesidedness of the mantle 
convection is in good agreement with the evidence given for these currents 
by the crustal deformations as shown by the data at the earth’s surface. 


1) Herskanen, W. A. and F. A. Ventne Merrnesz, The Earth and its Gravity 
Field, chapter 11, McGraw-Hill, New York—London (1958). Ventne Mernesz, F. A., 
The results of the development of the earth’s topography in spherical harmonics 
up to the 3lst order; provisional conclusions, Proc. Kon. Ned. Akad. v. Wetensch., 
Ser. B62, pp. 115-136 (1959); Convection-currents in the earth, Proc. Kon. Ned. 
Akad. v. Wetensch. Ser. B50, pp. 237-245 (1947); Gravity Expeditions at Sea, 
Vol. LV, Publ. Neth. Geodetic Commission, Delft (1948) 
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In the second place the current can only make a half turn. Because 
of the extremely small value of the temperature conduction by the mantle 
rocks, we can practically neglect the change of temperatures during such 
a half-turn current which lasts some 50-100 million years. We may, 
therefore, conclude that the temperature is carried along by the current 
and also the corresponding density deviations. The downward temperature 
gradient, brought about by the earth’s cooling— mostly by radiation at 
the earth’s surface —causes an unstable mass distribution in the mantle; 
the denser masses are situated in the top layer. In this statement we have 
neglected the density transition layer between 500 and 900 km depth 
which we shall presently discuss and of which we shall show that it does 
not affect our reasoning. 

The instability in the mantle by the earth’s cooling is stabilized by 
the elastic limit; an infinitely small disturbance can not bring it into 
action. If, however, some secondary temperature effect brings about a 
horizontal temperature gradient of sufficient size for causing a horizontal 
pressure gradient which can overcome the elastic limit, a convection- 
current must set in, gaining velocity until a quarter turn is accomplished. 
At that time the low temperature surface layer has moved into the sub- 
siding column of the current, and the high temperature lower mantle layer 
in the rising column; the two columns then have a maximum density 
difference and the current velocity has its greatest value. 

During the next quarter turn the density difference diminishes and 
disappears when a half turn has been made. The current must then come 
to a stop, or perhaps even somewhat before when the elastic limit can no 
longer be overcome. The mantle is then dynamically stable, and so we 
see that the convection-current must indeed have a half-turn character. 
The thermal conditions now are, however, far from stable ; the low tempera- 
ture matter is in contact with the core where the temperature is high, and 
the high temperature matter is in top where thus the temperature radiation 
towards the outside must be stimulated. During the period of rest, in which 
no mantle currents occur and, therefore, no great tectonic deformations 
of the crust, the temperature radiation and conduction must be especially 
active and this must on the long run bring back the original unstable 
density conditions in the mantle with low temperature on top and high 
temperature below. A new period of mantle convection with strong 
crustal deformations and consequently orogenesis must start. 

So we see that our mantle convection hypothesis, as already in 1939 
Grigas has pointed out, can explain the episodic occurring of orogenic 
periods, separated by long periods of rest. The periods of tectonic activity 
last about 50-100 million years, and this points to current velocities of 
a few inches per year. This checks with crustal movements during earth- 
quakes, which usually occur at intervals of a few tens of years. In this 
connection we may mention a relative shear movement during the earth- 
quake of 1892 in Sumatra of about 2.20 meters, and of 1906 in California 


28 


of about 4.50 meters. In both areas we know that these relative 
movements continued to occur in the same direction during long 
periods. 

The periods of rest lasted during several hundreds of million years and 
there is evidence that these periods have not been equal during the 
geological history of the earth. This last particularity is in harmony with 
the fact that the setting into motion of the mantle convection-currents 
depends on secondary processes which may, therefore, take place at 
arbitrary times. It is preferable, therefore, to use the term “episodic” in 
stead of ‘‘periodic”’ for the occurring of periods of tectonic activity during 
the earth’s history. 

Before coming to our main subject we have still to make a short mention 
of the problem whether the presence of the transition-layer at a depth 
from 500 to 900 km, where the density after reduction to surface tempera- 
tures and pressures increases downwards from 3.3 to 4.0, is compatible 
with our hypothesis of convection-currents through the whole mantle. 
This is the case if we follow an assumption already made in 1936 by 
BERNAL that the density transition is mainly caused by the change at 
greater pressure of the olivine into a spinel phase. A study of this problem 
was made by J. L. MetyErtne and C. J. M. Rooymans of the “Philips 
Laboratories” at Eindhoven!), and they arrived at a satisfactory con- 
clusion which can explain the great density difference and at the same 
time leaves open the possibility of convection-currents breaking through 
this transition layer; the driving power of such currents is in these con- 
ditions even considerably increased 2). According to MrtsprING and 
RooyMANS “the composition of the mantle probably deviates appreciably 
from the overall formula (Metal)gSiO4’’; enstatite is likely to be also 
present, and besides Fe and Mg several other metal ions. These sup- 
positions are entirely acceptable. 

We shall now make an attempt to attack our main problem concerning 
the question why the tectonic development on the east and the west side 
of the Pacific differs so much. For doing so we must realize that the 
presence of a system of convection-currents in the mantle must have two 
consequences. The drag forces exerted by these currents on the rigid crust, 
floating on the mantle, must in the first place bring about strong de- 
formations of the crust, which we shall presently examine. But in the 
second place we may expect that all these forces together exert a moment 
on the rigid crust as a whole, which tends to move the crust round the 
earth. In this connection the plastic properties of the mantle and the one- 


1) Meierrine, J. L. and C, L. M. Rooymans, On the olivine-spinel transition 
in the earth’s mantle, Proc. Kon. Ned. Akad. v. Wetensch., Ser. B 61, pp. 333-344 
(1958). 

ef) HEISKANEN, W. A. and F. A. Venting Metnesz, The Earth and its Gravity 
Field, § 11-2, pp. 406, 407; McGraw-Hill, New York, London (1958) 
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sided character of its currents is important; we can prove that in case 
the mantle consisted of a viscous Newtonian fluid, the moment of the 
drag-forces exerted on the crust would be zero. The crustal movement 
round the earth must combine with relative shifts of parts of the crust 
and with the crustal deformation which causes these shifts. The crustal 
deformation is often horizontal compression; in connection with our 
problem we shall confine ourselves to that case. Where the crustal 
deformation is caused by stress release we find graben or related phenomena 
with which we shall not deal in this paper. 

Horizontal uniaxial compressional stress of a size surpassing the elastic 
limit causes plastic downbuckling of the crust 1) leading to the development 
of a geosyncline. During the downbuckling stage of this development 
these belts are marked by belts of strong negative gravity anomalies, as 
have been found in the central parts of those island-ares which have been 
gravimetrically surveyed, as e.g. the Indonesian and Caribbean archipel- 
agoes 2). 

The downbuckling occurs in belts enclosing angles of about 55° with 
the direction of the uniaxial stress; this last direction may be considered 
to be identical with the direction of the subcrustal current causing the 
drag on the crust. On the west side of the Pacific these currents flow out 
from under the Asiatic continent, below which we may, because of the 
contents of radio-active matter of the continental sial, suppose a higher 
temperature than under the ocean and, therefore, the location of the 
rising columns of the convection-currents. 

It appears likely that, if no special conditions occur, the maximum 
uniaxial stress in the crust would be situated near the border of the 
continent, where usually the rigid crust is weakest; normally sediments 
are present there which can not have the same strength as the granitic 
and basaltic matter of the crust. The tectonic are is, however, much 
further away from the continent. Perhaps an explanation of this location 
might be found by assuming a gradual inward shift of the crust, with 
which the continental coast forms a whole, by a rotation of the crust 
round the earth. The geosyncline development, being caused by the drag 
exerted by the mantle current, may be supposed to remain more or less 
at its place. 


1) Ventne Mervesz, F. A., Plastic Buckling of the Earth’s Crust: The origin of 
Geosynclines; ‘““The Crust of the Earth”, Geol. Soc. 0. Amer. Spec. Paper 62, 
pp. 319-330 (1955). . 

2) Veninc Mernesz, F. A., Indonesian Archipelago: A geophysical study; Bull. 
Geol. Soc. 0. Amer., 65, pp. 143-164 (1954). Ventne Metnzsz, F. A. and W. A. 
HEISKANEN, The Earth and its Gravity Field, chapters 10A and 10C, McGraw-Hill, 
New York and London (1958). Venrnc Metnesz, F. A., Gravity Expeditions at Sea, 
Vol. IV, chapter II, publ. Neth. Geodetic Comm., Delftsche Uitgevers Mij., Delft 
(1948). Hess, H., Recent advances in interpretation of gravity anomalies and island- 
are structure. Adv. Report Comm. on continental and oceanic structure, (1939). 
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A similar phenomenon may be supposed to have occurred for all the 
island-ares east of Asia, the Marianas-Bonin arc, the Japanese arc, the 
Kuril arc, the Aleutian are, and also the Tonga-New Zealand are with 
other arcs east of the Australian-New Guinean continent. The crustal 
shift for all these arcs may more or less have had the same direction, 
viz. about N60°W, and the amount of the shift about the same, viz. 
about 15°. 

A certain support for this tentative hypothesis may perhaps be found 
in the fact that for several of these arcs other arcs may be found between 
the present are and the continent, which for some of them are known to 
be older and for the others may be surmised to be so. For the Indonesian 
are we may mention the Malay Peninsula~ NW. Borneo— Palawan — Luzon 
arc, the Sumatra ~ Banka — Bilitung —- Central Borneo ~ Sulu Islands — 
Philippine islands arc, and the Sumatra—Java-~Celebes -Sangi islands — 
Mindanao arc; for the Marianas - Bonin are at least two more ares to the 
west, beside the Nansei Shoto arc; for the Tonga - New Zealand are the 
New Zealand Ridge. The occurring of these intermediate ares between the 
present island-ares and the Asiatic continent is at least suggestive for the 
idea that in older periods the island-arcs east of Asia were situated at 
smaller distances from the continent. 

We shall now examine the east-coast of the Pacific. If the crustal shift, 
assumed above, has been world-wide, as it must be considered likely if 
we try to explain it by a rotation of the whole crust around the earth, 
the continents of North and South America must have covered the area 
where the mantle currents, rising under these continents and flowing off 
towards the Pacific, cause the greatest compressive stress. We may, 
therefore, expect that the compressive strain of the crust has occurred 
in the weakest belt of the crust in that area, viz. in the belt along the 
coast, without following the angles of 55° with the direction of the mantle- 
currents; we can understand that these angles can only be realized where 
the area under compression shows no appreciable differences of strength. 
We can thus account for the fact that on this side of the Pacific the 
mountain ranges follow the coast and that over large distances we find 
ocean troughs at the foot of the continental slope. In this way the great 
difference of the tectonic formations on both sides of the Pacific finds 
a natural explanation. 


Our hypothesis can also solve another problem which could hitherto 
not be satisfactorily answered, viz. the question what decides the shape 
of the island-ares. In the papers already referred to the writer has shown 
that each normal are consists of four parts, the two central parts making 
angles of 55°, to both sides, with the direction of the subcrustal current, 
where this current causes the crust to buckle downwards, and the straight 
parts on the two wings enclosing angles of 25° to 30° with this direction, 
where the main phenomenon is crustal shear though also some down- 
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buckling and some overriding of one side over the other is present. The 
lengths of these different parts, which decide the general shape of the 
island-arc, remained an open question. Our hypothesis solves this problem. 
The central downbuckling parts occur where the drag on the crust by 
the mantle current causes the maximum compression in the crust. The 
lengths of the wings depend on the amount of the crustal shift which 
withdrew the crust in the above mentioned direction of about N60°W. 
The general shape of the island-are is thus determined but details 
obviously depend on local circumstances. 
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A crystallographic investigation of the stable and unstable forms of some 
5-alkoxy-o.quinone-2-oximes and related compounds 
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(Communicated by Prof. E. Havinca at the meeting of November 28, 1959) 


Introduction 

In the course of the investigations performed in this laboratory on 
the tautomerism of quinone-oximes (nitrosophenols) [1] a study was 
made [2] of the two forms in which 5-alkoxy-o.quinone-2-oximes occur 
[3,4]. The absorption spectra give indications that in solution these 
compounds exist as tautomeric equilibrium mixtures, one of the isomers 
possessing a nitroso structure [2, 3,4]. This form which we call the 
f-isomer, predominates in non-polar media and can be obtained in a 
crystalline form from solutions in benzene, petrolether, carbontetrachloride, 
etc. The f-isomer transforms spontaneously into the x-form; the latter 
may be obtained directly by crystallisation from polar solvents such as 
alcohol. The «-form lacks the absorption at 600-800 mu, which indicates 
the presence of a nitrosogroup and which is found with the #-isomer 
(at 614-640 my). 

Although the reflection spectra and the infrared absorption spectra 
seemed to confirm the view that the «-isomer should be a quinone-oxime 
and the f-isomer a nitrosophenol derivative, a crystallographic analysis 
was deemed essential to give independent direct evidence on the tautomeric 
relationship of the «- and the f-forms and to obtain data on the hydrogen 
bond formation between the molecules, the syn- or anti-configuration, 
and the conformation of the side chain. To this end one of us extended 
the number of substances available for investigation by synthesising the 
5-(2'-bromoethoxy)- and the 5-(2’-chloroethoxy)-o-quinone-2-oxime. These 
compounds again were found to occur in two forms: stable orange-red 
needles and unstable green lath-shaped crystals. 

In this paper the crystallographic analyses of the compounds and the 
conclusions on the steric structures of the isomers will be reported 1). 


The stable forms 
By slow evaporation of alcoholic solutions of «-5-(2’-R-ethoxy)-o.quinone 


2-oxime (R= Br, Cl) orange-red needles are obtained. The data summarized 


') Details of the syntheses and a discussion of the absorption and the reflection 
spectra can be found Le. 2). 
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Spacegroup, dimensions (Angstrém units) of the unit cell and densities of 
x-5-(2’-R-ethoxy)-o.quinone-2-oximes at room temperature 


R a b c B n d 


‘obs doale spacegroup 
Br 4.34 20.2 HA AY) 105 4 P 21/e 
Cl 4.13 19.62 1S, 105.4 4 1.62 1.54 P 2i1/e 


in table 1 show that the bromoethoxy- and the chloroethoxy-compound 
are presumably isomorphous which assumption was confirmed by sub- 
sequent analysis of the X-ray diffraction data. The a-axis was taken in 
the direction of the needle axis. 

With the aid of well-known isomorphous replacement procedures the 
well-resolved projection along [100] of the structure of the «-chloroethoxy- 
compound was studied. From Okl reflection data the normal projection 
oo (y 2) was calculated (fig. 1) and from 1k] reflection data the generalized 
projection 0; (yz). The stage of refinement with a reliability index of 
16 % for 0kl structure factors and of 23 °% for 1k] structure factors does 


-Vob -Yy,b 0) Yb Yb 


Fig. 1. The electron-density 0(y 2) of «-5-(2’-chloroethoxy)-o.quinone-2-oxime, 
projected along [100]. The contour lines are drawn at .2, .4,... eA-2, except those 
for the chlorine atoms which are drawn at .2, .5, 1.0... ceA-2, The .2 eA-2 line is dashed. 


not permit a definite decision in favour of an o.quinone-monoxime 
structure or of a nitrosophenol tautomer. However, the quinone-monoxime 
structure has been established beyond doubt for x-2-chloro-5-methyl- 
quinone-4-oxime [1, 5] and the same is true for 3-chloro-quinone-4-oxime 
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[1] and for the first representative of the homologous series, «-5- 
methoxy-o.quinone-2-oxime [6]. 

In all these compounds hydrogen bonds which connect the carbonyl 
group of one molecule with the oxime group of the next molecule act as 
links in infinite chains in the crystal. «-5-(2’-chloro-ethoxy)-o.quinone- 
2-oxime exhibits the same kind of hydrogen bonds. Its absorption spectrum 
does not show the typical nitroso absorption band [2]. It is therefore 
probable that the «-form has an 0.quinone-monoxime structure with the 
NOH group in anti-configuration with respect to the CO group. Of special 
interest is the Cl-C2H4-O group. The chlorine and the oxygen atoms of 
this group are present in a gauche conformation. 


The unstable forms 


The f-modifications are obtained by slow cooling of saturated solutions 
in non-polar solvents. The chloroethoxy-compound is very unstable; we 
met several difficulties in obtaining suitable crystals and subsequent 
X-ray diffraction pictures. A comparison of the unit cell dimensions 
(table 2) and of the intensities of the diffraction spots proved the 
B-bromoethoxy- and £-chloroethoxy- as well as the £-n.propoxy-compound 
to be isomorphous. Until now the analysis has been focused on the 
B-chloro compound, since the measurements of the cell dimensions of the 
f-n.propoxy compound have been only recently determined. 


TABLE 2 


Spacegroup, dimensions (Angstroém units) of the unit cell and densities of 
B-5-(2’-R-ethoxy)-o.quinone-2-oximes (nitrosophenols) 


R a b c B nm Aone degie  Spacegroup temperature 
Br 6.76 12.73 10.86 90 4 41.77 # £1.75 P 2i/a 20°C 
Cl 6:69 12.28) 10,70 4 30) —4)- 9148) 91.64 P 21/a —140°C 
CH, 6.81 12,63 10:61 “S065 14) > 1:20 126 P A/a —140° C 


The isomorphous replacement procedure was applied and the normal 
projection go (yz) and the generalized projection 0; (y z) were calculated. 
The resolution of the projection along [100] is quite good (fig. 2), with 
the exception of the carbon atom belonging to the carbonyl group. This 
atom overlaps itself by the action of a center of symmetry. 

The molecules are found in layers perpendicular to the needle axis 
at heights x=} and a=. Therefore, we calculated the “mixed” projection 
00 (y 2) +Si (yz) based upon structure factors F(0kl) and F(1kl), respect- 
ively 1). S denotes the imaginary part of the generalized projection [7]. 
In this way (fig. 3) the above-mentioned carbon atom became resolved. 
The result of the analysis can be summarized as follows: 

1. there are no intermolecular bonds; 


‘) We only succeeded in measuring intensities of 0kl and 1kl reflections of the 
B-chloroethoxy-compound. 
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the NO(H) group is syn with respect to the CO(H) group; 
an intramolecular hydrogen bond is probable; 

the R-CyH4-O group is in gauche conformation: 

except for the R-C:H4-O group the molecule is planar ; 


a RS ee oa 


the crystals of the f-form are layer structures. 


Y 
ab 


-4b 
— \gcsinB 0) 1/2 ¢ sinB Sig ¢ sinp 
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Fig. 2. The imaginary part Si(yz) of the generalized projection oi(yz) of 
B-5-(2’ chloroethoxy)-o.quinone-2-oxime (nitrosophenol), projected along [100]. The 
contour lines are drawn at arbitrary intervals. The contour lines belonging to the 
chlorine atoms are drawn on a scale 2.5 times as large as those of the light atoms. 
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Fig. 3. The “mixed” projection g,(y z) + Si(y z) of B-5-(2’ chloroethoxy)-o.quinone- 

2-oxime (nitrosophenol) projected along [100]. The contour lines are drawn at 

arbitrary intervals .The contour lines belonging to the chlorine atom are drawn on 

a scale 2.5 times as large as those of the light atoms. The lengths of the bond 
distances are indicated. 
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The refinement of the structure stopped at the stage R=16 % for Okl 
reflections and at the stage R=20 % for 1kl reflections. The approximate 
dimensions of the molecule are given in fig. 3. 

We estimate the distances to be accurate to about 0.08 A. The ring 
system of six atoms is a nearly regular hexagon with mean C-C distances 
of 1.40 A. The C—O distance (1.28 A) is quite large for a carbonyl group, 
but is consistent with the presence of a COH group in a mesomeric system 
in which the nitrosophenolic limiting structures are predominant. These 
bond distances do not establish the nitrosophenolic character of the 
6 form beyond doubt; but the optical properties of its crystals support 
the impression that the 6 form is mainly a nitrosophenolic structure. Its 
powder spectra show the nitroso absorption band at co 610 my. This 
absorption is due to a n—z™* transition [8, 9, 10] and is perpendicular to 
the plane of the molecule. 

As stated before the molecules are equidistant in planes a=} and 
x=# and are perpendicular to the a axis. A very marked pleochroism [3] 
is found by inspecting the crystals in polarized light. The crystals are 
orange if the needle axis is perpendicular to the electric vector of the 
polarized light, but grass-green if the needle axis is parallel to the 
electric vector. 

This case closely resembles the example of pleochrosim occurring in 
crystals of m-nitronitrosobenzene given by Nakamoto and Suzuki [11]. 
These authors could give direct evidence that the n—z* absorption is 
perpendicular to the plane of the molecule. 

Crystals of monomeric nitroso compounds e.g. p.iodonitrosobenzene 
[12, 13] that have the typical absorption band at 700 mu show pleochroism. 
But in the latter case the pleochroism is a not very pronounced shift 
between lighter and darker shades of green, because the molecules are 
not packed in layers. 

We intend to continue the study of the crystal structures of «-5-chloro 
(2’-ethoxy)-o.quinone-2-oxime and of B-n.propoxy-o.quinone-2-oxime. The 
authors are indebted to Dr. J. C. Schoone, State University of Utrecht, 
who made structure factor calculations on the Z.E.B.R.A. They 
wish to express their gratitude to Prof. Dr. E. Havinga for his stimulating 
interest. 
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PALEONTOLOGY 


MIOGYPSINA IN NORTHWESTERN GERMANY 
BY 


C. W. DROOGER 


(Communicated by Prof. G. H. R. von KoENIGSWALD at the meeting of Nov. 28, 1959) 


ABSTRACT 


Two species of larger Foraminifera have been found in the classical deposits of 
Astrup and Doberg. Miogypsina septentrionalis n.sp. clearly belongs to the M. 
complanata group, which is of Late Oligocene age in the Mediterranean area. The 
few accompanying Lepidocyclinids are closest to Lepidocyclina morgani of the same 
southern area. Both species enable us to make a direct comparison with the sedi- 
mentary sequence of the Aquitaine basin, proving that the classical Chattian 
deposits of the North Sea basin are older than those of the typical Aquitanian near 
Bordeaux. 


INTRODUCTION 


In 1956 samples were collected from the Late Oligocene ‘“Meeressande”’ 
of northwestern Germany: near Kassel (Ahnetal), Osnabriick (Astrup) 
and Biinde (Doberg). In some of our samples small numbers of larger 
Foraminifera were encountered. Their extreme scarceness and rather bad 
preservation may account for the fact that they have not been discovered 
earlier. In the papers of GoLtpruss, Reuss and others we found no clear 
indications that these authors had classified these Foraminifera under 
the Bryozoa. 

The localities Kassel, Astrup and Doberg are famous for their 
fossiliferous deposits, the faunas of which have been studied by many 
authors since the first half of the 19th century. With the sediments near 
Kassel as the type, Fucus (1894) introduced the Chattian, the now 
internationally adopted member of the general time scale for the Late 
Oligocene. Foraminifera from these localities were described as early as 
1838 (RopmpR). Several additional notes and papers were published in 
later years, and Batsns (1958) has recently given a thorough revision 
of many of the species. 

Outcrops of such Upper Oligocene calcareous microbreccias with many 
molluses, echinoids, bryozoans, ete., occur at a few places only, widely 
scattered in northwestern Germany. These outcrops are probably remnants 
of a formerly much greater area of extension of shallow water sediments. 
In subsurface sections the corresponding beds are less calcareous and less 
fossiliferous (ScHAD, 1947). In some more recently investigated sections 
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in the Lower Rhine area they were found to contain rich microfaunas 
(InpAns, 1958; ELLERMANN, 1958). Very similar associations from cor- 
responding beds have been described from the Netherlands’ Peel region 
(TEN Dam and REINHOLD. 1942) and from northern Belgium (BatTJEs, 
1958). These deposits have been given several names: lately they were 
united under the name Elphidiwm-Almaena beds, because of their 
characteristic species Elphidium subnodosum (RoEMER) and Almaena 
osnabrugensis (ROEMER) (DROOGER and Barszs, 1959). 


LOCALITY DETAILS 


The samples taken during the trip of 1956 were derived from: 

1. the Doberg near Biinde. In a small syncline some 70 metres of 
more or less indurated, calcareous microbreccias have been preserved; at 
least this measure of thickness was visible at the time that Hupacu 
established the detailed section (1922), which was published many years 
later, in 1957. At the moment the lower part is no longer visible. HuBacH 
distinguished two parts in the section, which were published in 1957 
under the names: Lower and Upper Doberg beds. This nomenclature is 
considered preferable to that of HuBAcH’s manuscript of 1922 (Eochatti- 
cum and Neochatticum), which unfortunately was adopted by some later 
authors, and which may cause an unnecessary confusion of time and 
rock units. 

The Pectinidae of the Doberg appeared to be very suitable for strati- 
graphic subdivision, allowing the distinction of several biozones, which 
could be recognized at other localities as well. Already Husacu distin- 
guished a phylogenetic series in this group of molluscs. Later authors 
(GorGES, 1951-1957; ANDERSON, 1958) were of different opinions on the 
determinations, but the basis for the subdivision remained. Unfortunately, 
these species of the Pectinidae seem to be restricted to the North Sea 
basin. Within this basin the Doberg appears to offer the most complete 
exposed section of Upper Oligocene rocks. 

Some ten samples were taken throughout the visible Doberg series, 
most of them from the Upper beds. Only one sample, from friable, 
yellowish, fairly homogeneous, fine debris with much glauconite (DB 3), 
was found to contain Miogypsinids. There is an accompanying rich 
microfauna, with large-size species, including Elphidiwm subnodosum 
(RoEMER), Planorbulina difformis RopMER and many Lagenids, such as 
Frondicularia oblonga (RormmER), Vaginulinopsis gladius (PHILIPPI), 
Sigmomorphina regularis (RomMER) and Guttulina problema D’ORBIGNY. 
Macrofossils are scarce; there are numerous bryozoans, very small 
echinoids, and some immature Pectinids. The sample was taken in the 
northern wing of the syncline from one of the lowermost visible beds, 
at the level of the road and north of it, and in front of the open air theatre. 
Dr. ANDERSON kindly informed us that the sample had been taken from 
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the higher part of the Lower Doberg beds, in the zone of Chlamys 


hausmanni (GOLDFUSS). 

2. Astrup, the abandoned marlpit, about 2 km NNE of Belm near 
Osnabriick (see BarsEs, 1958, p. 69). There is a discontinuous exposure 
of up to 5 metres height in the roughly N-S striking western wall of the 
former excavation. The wall shows several normal faults with unknown 
amounts of displacement. As a consequence the exact spacing of the 
samples from different blocks is unknown. Three of the samples were 
found to contain larger Foraminifera. 

BA 1 was taken in the most southern block from coarse greenish sand 
with much calcareous debris and glauconite. About 50 em lower and 
150 cm higher sandy gravels begin, containing pebbles of up to 10 cm 
diameter. Some 20 Miogypsinids were encountered. 

BA 3 was taken at the level of the bottom of the exposure, some 
10 metres further north in the next block, which has a downthrow relative 
to the more southern block of unknown amount, but it is most likely 
that BA 3 is from a higher level than BA 1. The sample was derived 
from whitish, fairly homogeneous, fine, caleareous microbreccia with an 
admixture of but rare other components, mainly glauconite. Some 2} m 
higher there is a somewhat indurated layer of about 15 em thickness 
with numerous casts of molluscs. Again some 20 Miogypsinids, and a 
few Lepidocyclinids were found. 

BA 8 came from the southern side of the hollow road which passes 
along the most northern tip of the former pit’s wall. There was no 
continuous exposure between BA 3 and BA 8. The sample was taken 
some 1$ metres above the road’s level from coarse, ill-sorted, greenish, 
calcareous sand with molluscan fragments and other debris, glauconite 
and small, usually flattened, gravel of variable nature. Higher up, the 
sand and gravel become coarser with elements of up to 6 cm diameter. 
Some fragmentary Lepidocyclinae are the only larger Foraminifera. 

All our other samples from Astrup were devoid of larger Foraminifera. 
Especially BA 3 contains many bryozoans and a fairly rich microfauna 
of about the same composition as DB 3 of the Doberg. According to 
Dr. ANDERSON (personal communication) the layers exposed in this wall 
are comparable to the Upper Doberg beds, the zone of Pecten hofmanni 
GOLDFUSS. 

3. Kassel. By courtesy of the Hessischen Landesamt fiir Boden- 
forschung (Wiesbaden), a fresh exposure had been made in the Ahnetal, 
below the Brandkopf. At the time (1956) this exposure did not show the 
underlying Septaria clay. Some 44 m of mainly yellowish to greenish, 
calcareous sands were visible. All our seven samples contained macro- 
fossils, but no larger Foraminifera. The foraminiferal associations show 
much resemblance to those of the previous localities, but the big Lagenid 
specimens are mainly accompanied by numerous Miliolidae, Bryozoan 
debris appeared to be much less frequent. 


4] 
SYSTEMATIC DESCRIPTION 
Altogether over 50 Miogypsina specimens were available. Only some 


eight fragments of Lepidocyclina were met with, most of them in a very 
bad state of preservation. 


Family Miogypsinidae Tan Sry Hox, 1936 
Genus Miogypsina Sacco, 1893 
Subgenus Miogypsina Sacco, 1893 
Miogypsina septentrionalis n.sp. 


Pl. 1, fig. la, b, 2; pl. 2, fig. 1-8 


The species was given the name septentrionalis because of its far 
northern occurrences, which are exceptional within the genus. 

The holotype (coll. no. $ 12165, of the paleontological collections of the 
Mineralogisch-Geologisch Instituut, State University of Utrecht) was 
selected from sample DB 3, from the upper part of the Lower Doberg 
beds (Chattian) of the Doberg near Biinde, Westfalen. Other occurrences 
of the species are in samples BA 1 and BA 3 from the Upper Oligocene 
(Chattian) of Astrup near Osnabriick, Niedersachsen. The individuals 
of the Doberg are slightly less well-preserved than those from Astrup, 
but because of the clearer stratigraphic position of the former locality 
sample DB 3 was taken as the type sample. 

Diagnosis. The species is characterized by the high average number of 
nepionic chambers (about 16), the low average number of operculinid 
spiral chambers (about 3 to 4), and the development of lateral chambers, 
of which there is usually not more than one layer. It is different from the 
morphologically nearly identical, but probably not directly related 
M. thalmanni DrooGER mainly in the number of operculinid chambers 
(average 8 to 9 in M. thalmanni), and from some assemblages of the 
closely related M. complanata SCHLUMBERGER 8.1. mainly by the presence 
of lateral chambers. 

Description. Nearly all individuals of DB 3 are damaged, especially 
along the frontal edge. The apical border is usually better preserved and 
more or less entire. The general shape of the test is roughly circular to 
fan-shaped, but it is often obscured by the considerable damage. The 
diameter of the test ranges from 1.0 to 1.8 mm, the thickness from 0.4 
to 0.6 mm. Greatest thicknesses occur in the broadly rounded apical 
portion over the early spiral part. The individuals are often unequally 
biconvex. 

The surface shows flat pustules of 50 to 125 w in diameter. Usually 
they are of the fairly uniform size of about 100 p. The appearance of 
several specimens gives the impression that they had been rolled off and 
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slightly corroded, so that the pustules have become indistinct. In most 
individuals the surface shows distinct holes at irregular places in between 
the pustules as if the roofs of some lateral chambers had been removed 
or had not been formed. 

During the grinding the median layer appeared to be variously curved, 
which made it often impossible to obtain good sections of the entire 
median layer. Altogether ten satisfactory sections can be used for the 
analysis of the internal characteristics. 

The protoconch, of fair size, was found to be more or less circular in 
section. The deuteroconch, measured at right angles to the line connecting 
the centres of both embryonic chambers, is of about equal diameter. 
The single nepionic spiral consists of one to one and three quarters of a 
coil. On the average, increase in size of the spiral chambers continues 
up to the tenth chamber (10.5 + 1.2; range 18 to 4), before decreasing 
again, but the individual observations are often far from accurate, 
differences in size of successive chambers commonly being hard to observe. 
A similar difficulty of observation concerns the point of beginning of 
orbitoidal growth, the first equatorial chambers being very small and 
situated towards the apical border. Operculinid growth continues during, 
on an average, only four nepionic chambers (4 + 0.5; range 6 to 1). In 
several specimens these first equatorial chambers are intercalary chambers, 
i.e. they are covered by the final chambers of the main nepionic spiral. 
Some observations (especially in BA 3) showed that the nepionic spiral 
is slightly trochoid. 


Internal features of Miogypsina septentrionalis 


Sample N be y° diam. I (u) 


nr. R M + oy R M + oy R M + oy Mal Mr 


22 — 11/16.0 + 1.0 |-270 — -140|-215 + 14.4] 140— 180} 160+ 4.8] 0.98 
BA 1 | 3/21— 15/18 360 — —210|-—280 140 — 180} 160 
BA 3 |10/20— 13) 15.9 + 0.7 |-260 — -120| -190 + 17.0) 130 — 210 164 + 8.1 0.92 


N = number of specimens. 
X = number of spiral chambers following the deuteroconch. 


y = angle between the apical-frontal line and the line connecting the centres of 
both embryonic chambers. 
I = protoconch, II = deuteroconch. 


R = range, M = mean. 


Equatorial chambers are relatively few in number; up to five rows of 
them were observed. Their shape in horizontal section is arcuate to 
ogival; maximal observed dimensions are 175 x 225 Hu. Occasionally they 
show distinct basal stolons. 

Transverse sections as well as broken individuals show that the side 
walls contain distinct lateral chambers. There is usually not more than 
a single layer on either side, but as many as three lateral chambers were 
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occasionally observed on top of one another. The layer is commonly 
discontinuous, mainly because of the irregular shape of many of the 
chambers in cross section. Some give the impression of being vacuoles, 
originated in the septa between the spiral or other median chambers, 
but several observations confirm that the majority are true lateral 
chambers, often connected with the median chambers through stoloniferous 
openings. 

The individuals of the Astrup samples are hardly different from those 
of the Doberg. They are somewhat better preserved, and a slightly greater 
diameter was measured (up to 2.0 mm). Correspondingly the maximal 
number of rows of equatorial chambers is greater, up to six. The nearly 
identical values for the internal features are shown in the table. For 
BA 3 the average number of operculinid chambers is 3.3 + 0.5 (range 
7 to 3), of increase 11 + 0.8 (range 16 to 8). Again, lateral chambers 
appear in small numbers on either side of the median layer. 

Remarks. The presence of distinct lateral chambers, though few in 
number, place our assemblages in the subgenus Miogypsina. Among the 
species of this subgenus our German assemblages are exceptional on 
account of the the high average numbers of nepionic chambers. Only 
Miogypsina thalmanni DROOGER (1952, pp. 15, 47, pl. 1, fig. 1-5, pl. 2, 
fig. 1-5, pl. 3, fig. la, b; table 1) from the Oligocene Tabera formation 
of the Dominican Republic, has similarly high My values (15.2), but 
the average number of operculinid chambers is higher (8.5; R=13-5; 
iV =18) and the protoconch size much smaller (My;=90 vw; R=60-115 yp). 
The increase of size of the nepionic chambers is about equal (= 10.7; 
P= -7): 

The groups of populations, those from Germany and those from the 
Dominican Republic, are both very close to the origin of Miogypsina 
s.str. in their own region, but because of their wide geographic separation 
there seems to be little doubt that they are different stocks, which evolved 
independently from Miogypsinoides. Since the European assemblages can 
morphologically be clearly distinguished on features of the early stages, 
a separate specific name has a sound basis. 

There is also great morphologic resemblance to several populations of 
M. (Miogypsinoides) complanata SCHLUMBERGER s.l. For correlation 
purposes the assemblage from the Ferme Christus near Dax, SW France, 
is very important. Apart from the absence of distinct lateral chambers 
in the French individuals, there is a fair resemblance. The features of the 
early stages are: N=15; Ry = 21-9; My =13.9+ 0.74; R,= —270°— — 80°; 
M,= —175° + 16.1; Rr=110-225 w; Mr=153 w + 9.4; My/My=1.03 
(DRooGER, KAASSCHIETER and Kny, 1955, table 1); operculinid chambers: 
R=7-2; M=4.3; chamber increase: R=12-4; M=8.2. In nepionic 
features the assemblage from Christus seems even to be slightly more 
highly developed, but differences are very small. Another assemblage 
very close to the German ones is that from the Bric del Duca (sample 2) 
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in northern Italy (DRouGgER, 1954, table 1), but this assemblage is still 
more advanced in its nepionic features. 
There is very good reason to suppose that some stock of M. complanata 
3 very 2 : 
of average M, value migrated from southern Europe to the North ae 
basin, where it developed somewhat different features, by which it crosse 
the morphological boundary of the subgenus Miogypsina. 


PLATE 1] 


All figures x 30. 
Miogypsina (Miogypsina) septentrionalis n. sp. Upper Oligocene. 
la, b. holotype. a, side view; b, peripheral view. DB 3, Doberg. 2: side view. BA ay 
Astrup. 


Figs. la, 6, 2. 


Figs. 3, 4. Lepidocyclina (Nephrolepidina ) morgani Lemotné and R. Dovvirxe. 
Upper Oligocene. 3. side view. BA 3, Astrup. 4. peripheral view of an irregularly 
broken specimen. The left half of the figure shows a transverse section near the 
centre of the test, and the protruding glauconite filling of the embryonic chambers. 
The right half is much more peripheral. The median layer in between is obseured. 


BA 8, Astrup. 
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Family Lepidocyclinidae Scurrren, 1932 
Genus Lepidocyclina Gimprn, 1870 
Subgenus Nephrolepidina H. Dovyi1#, 1911 
Lepidocyclina morgani Lemorne and R. Dovuvim1t 


PL, fig. 3,4, pl. 2; fig. 9, 10 


Lepidocyclina morgani LemoinE and R. Dovvitrt, 1904, Soc. géol. France, 
Paléont., mem. 32, p. 17, pl. 1, fig. 12, pl 2, fig. 4, 12, pl. 3, was 2; DRooGER and 
Soctn, 1959, Micropaleont., vol. 5, p. 420. 


Lepidocyclina tournoueri, eee (not L. and R.D.), 1940, Schweiz. Pal. 
Abh., vol. 63, pp. 47-54, pl. 3, fig. 3-6, 15, pl. 5, fig. 3, 11, 12, 14; Droogmr and 
Macne (not L. and R.D.), 1959, Micropaleont., vol. 5, p. 279, textfig. 4. 


Only three fairly complete specimens and a couple of fragments were 
available, all from Astrup, samples BA 3 and BA 8 

The shape is regularly lenticular with observed diameter between 
1.7 and 2.1 mm, and a thickness of about 0.7 mm. The surface shows 
flat pustules all over with numerous irregular hollows in between. This 
surface possibly has been modified by later partial solution of lime and 
other damage. 

Horizontal sections of two specimens fairly well show the features of 
the median layer. The embryonic stage is nephrolepidine and of nearly 
identical size in both individuals. The protoconch has a larger diameter 
of about 350 yw, the deuteroconch of about 500 w, and both embryonic 
chambers together, measured through their centres, of again roughly 
500 yw. There are symmetric nepionic spirals along the protoconch. The 
deuteroconch shows three or four accessory auxiliary chambers along its 
circumference. The shape of the equatorial chambers varies from arcuate 
to ogival; in the external part these chambers are more or less arranged 
according to a pattern of intersecting curves. The observed dimensions 
near the border are 75 x 75 wp. 

Some fragments show that the lateral chambers are flattened, much 
larger than the equatorial chambers, and arranged in columns, those of 
adjoining columns more or less alternating. 

The observed features indicate that our specimens are closely related 
to the group of Nephrolepidina populations of the Mediterranean area, 
which belong to the L. morgani-tournouert group (see also VAN DER VLERK, 
1957). No sound conclusion can be drawn from our two sections, but the 
number of accessory auxiliary chambers, considered to be the primary 
diagnostic feature in this group, points to L. morgani. Assemblages from 
the Aquitaine basin, such as those from Saint-Etienne-d’Orthe, Saint- 
Géours and Abesse (unpublished observations), in which L. morgani 
is accompanied by M. complanata s.l., invariably show averages for the 
number of accessory auxiliary chambers between three and four. Also 


PLaTEe 2 


All figures x 30. 
Figs. 1-8. Miogypsina (Miogypsina) septentrionalis n. sp. Upper Oligocene. 
1-4, transverse sections showing the slight development of lateral chambers. Only 
section 4 passes through the outer part of the nepionic spiral, the others are more 
nearly cut along the apical-frontal line. 1 and 2 from BA 3, Astrup; 3 and 4 from 

DB 3, Doberg. : 
5-8. median sections; section 6 has not been drawn completely. 5 from BA iF 
Astrup; 6 and 8 from BA 3, Astrup; 7 from DB 8, Doberg. 
Figs. 9, 10. Lepidocyclina (Nephrolepidina) morgani LEMOINE and R. Dovvittez, 
Upper Oligocene. Median sections, of which 10 has been drawn only partly. 

BA 8, Astrup. 
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the ratios between the diameters of I and II are similar. Only the absolute 
values of the diameters of the embryonic chambers of our two German 
specimens should be extreme and high in the variation pattern of the 
French group. 

Again, some regional differentiation may have developed in our North 
Sea populations, but the main features make as yet the specific determi- 
nation as L. morgani the most likely one. 


CONCLUSIONS 


Possibilities of close comparison between the typical Chattian and 
Aquitanian sediments have been few. This resulted in the numerous 
discussions and very diverging opinions about their relative stratigraphic 
position. The North Sea basin and the Aquitaine basin are wide apart, 
without continuous good connections during Oligocene and Miocene 
time; moreover, sedimentary and climatic conditions were probably 
much different in both basins during most of this time. 

The fact that at least one wave of larger Foraminifera of two independent 
and fairly well-known lineages penetrated from the Mediterranean area 
into the North Sea basin, gives a sounder basis to correlations. 

According to the Pecten zoning of the northern area, the sediments 
with larger Foraminifera of Astrup and Doberg are slightly younger 
than the type deposits of the Chattian near Kassel. However, there is 
general agreement about their incorporation in this Late Oligocene stage. 
The type sediments of Kassel are comparable to the middle part of the 
Lower Doberg beds. The absence of larger Foraminifera in our samples 
from Kassel is probably due to environmental circumstances. The 
abundance of Miliolidae points to less open-marine conditions during 
sedimentation in the Kassel area. 

Miogypsina septentrionalis has been shown to be closely related to, and 
to be a local somewhat divergent stock of M. complanata of the Mediter- 
ranean. Because of the principle of nepionic acceleration there seems to 
be little doubt that it is also comparable in age. This seems to be confirmed 
by the scarce data on the accompanying Lepidocyclina morgant. 

Associations of M. complanata and L. morgani with the exclusion of 
other species of these groups are fairly common in the Mediterranean area. 
Most important for stratigraphic relations are those of the Aquitaine 
basin, where they are known from the region of Dax in the southwest, 
from classical localities such as Saint-Etienne-d’Orthe, Saint-Géours-de- 
Maremme (Escornebéou) and Abesse. 

There has been much discussion about the age determination of these 
deposits. Specialists of molluscs had varied opinions, even ranging from 
Tongrian to Helvetian. In later years a Late Oligocene age was thought 
most likely (Dacuty, 1948), but particularly among micropaleontologists 
the majority of the localities were considered to show sediments of 
Aquitanian age. This certainly happened under the influence of the 
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authors on larger Foraminifera (LEMOINE and R. Dovuvim1h, 1904; 
H. Dovviit, 1925), who distinguished several zones, mainly on the 
basis of relative proportions of, and presence or absence of, the represent- 
atives of the various genera and subgenera. As most localities involved 
did not contain representatives of Nummulites, they were believed to be 
younger than Oligocene (H. Douviii6, 1907). Abesse was placed in the 
Burdigalian, St. Etienne-d’Orthe in the Burdigalian or Aquitanian, and 
St. Géours in the Aquitanian or Latest Oligocene. The latter determination 
was made by H. Dovvit1h in 1925, who reported small Nummulites 
from this locality. We failed to find representatives of this genus at 
St. Géours, but no doubt the exposures have become much worse since 
most marlpits have been abandoned. The still visible part at St. Géours 
(Escornebéou) was again placed in the Aquitanian by Sian (1949) and 
other authors. 

It must be emphasized that close comparisons were never made with 
the faunas of the type deposits of the Aquitanian near Bordeaux. 
Investigation of the microfauna (DRooGER, KAASSCHTETER and Key, 
1955) has shown that the typical Aquitanian sediments contained more 
highly developed species of the Miogypsina lineage (M. gunteri CoLE 
and M. tani DroogerR) and that they were devoid of Lepidocyclina. 
Hence, it was concluded that all deposits of the southwestern Aquitaine 
basin with the associations of M. complanata, of various heights of 
development, and L. morgani (erroneously named L. touwrnoueri) were of 
pre-Aquitanian age. 

The discovery of nearly the same association in sediments of undoubtedly 
Chattian age of the North Sea basin fits in nicely with these more recent 
opinions about the deposits of the Aquitaine basin. It may now be 
safely decided that the Chattian and the Aquitanian are 
really successive time intervals without notable overlap 
and without a considerable gap in between, since M. complanata 
s.l. of the Chattian and M. gunteri of the Early Aquitanian are successive 
species of the main Miogypsina lineage in the Aquitaine basin. The less 
well-known lineage of European Nephrolepidina thus far confirms this 
opinion, as our few North Sea specimens point to a resemblance to the 
Chattian populations of L. morgani and not to those of the Burdigalian 
L. tournouert or to as yet unperfectly known populations in between. 

Among the smaller Foraminifera another peculiar species of both basins 
now gets a much more satisfactory stratigraphic distribution. Almaena 
osnabrugensis (ROEMER) with all its local variants and possible subspecies 
has been frequently reported from the sediments of the southwestern 
Aquitaine basin under discussion, such as those of Escornebéou, where 
it 1s very abundant (SiaaL, 1949, 1950; HorKkmr, 1952; Rey, 1958). 
According to our new correlation it would not range beyond the Chattian 
of the Aquitaine basin, which is exactly the condition in the North Sea 
basin, where it is also restricted to the Oligocene. 
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PEBYSIICS 


ON THE POWER TRANSFER BETWEEN PARAMAGNETIC SPINS 
AND CRYSTAL LATTICE. IIIa 


EXPERIMENTAL RESULTS 
BY 
B. BOLGER 


(Communicated by Prof. C. J. Gortser at the meeting of June 27, 1959) 


A. Results obtained by steady state saturation 
[3.1] Presentation of the data 

The steady state saturation measurements were all made at a frequency 
around 9400 MHz. Graphs have been made of P,,,/Pm vs P,,, a8 calculated 
from the measurements with formulas (2.15) to (2.17) or (2.24). The 
relation of P,.,/Pm to better known quantities is that 


Pest 1 
joy ome be ocers 


Pm is taken for one gram ion of magnetic ions. 7, the molar power transfer 
constant, is the reciprocal of the slope of the graphs 

n= dPext 

dP. ext/ is m) 

7” was expected to be a constant but mostly, however, this is found not 
to be realized. Therefore 7 at low values of P,,, is called ys, and for high 
values 7. These as well as the coordinates of the saturation graphs will 
then be given. Also are given the coordinates of the point at which 7 = 2ns 
in the column “coord. 7=27,”. When the range of P.., was not large 
enough to observe a definite curvature for a series no index will be 
assigned to the corresponding 7. To convert the values of 7 to relaxation 
times one has to use (1.89) with W=0 in case of internal spin equilibrium. 

Since the saturation graphs were expected to be straight lines instead 
of the curved ones found, different experimental methods (reflection and 
transmission) with different cavities have been tried out to exclude an 
instrumental error. Also the apparatus was tested by experiments on 
CuK tutton salt and CrK alum diluted with AIK alum. Apart from CuK 
tutton salt at the lowest power levels, nice straight lines were obtained 
for these salts up to a very high degree of saturation, indicating that no 
instrumental error was involved. 

Where possible a comparison is made with the results obtained by 
other experiments, especially with those obtained by the (non-resonant) 
relaxation method. The values of 7 (called 77,,,) derived from those results 
are calculated on the assumption that the spin system is in internal 
equilibrium, using (1.91). 

Our measurements have all been made on single crystals while the 


(3.01) 
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relaxation measurements were mostly performed on powders. As in 
general 7 is expected to be anisotropic, the main point of comparison will 
be the temperature dependence and the order of magnitude. 

For each salt studied the crystal structure and magnetic data are given. 
The crystal axes will be denoted by a, b and c, the magnetic ones by 
z, y and z. 


[3.2] Results on CrK alum 

Experiments were carried out on the trivalent Cr**+*+ion in single 
erystals of KCr(SO4)2 12H2O, some of the experiments being carried out 
on crystals magnetically diluted with the isomorphous AI salt. 


[3.2.1] The crystal structure of the alums has been investigated, 
in particular by Lipson and Brrvers [1] [2]. The crystals are cubic 
(space group Pa 3) with four molecules in a unit cell. More detailed analysis 
showed that one can distinguish at least three types of structures the 
«x, B and y structure. Which of them is the stable one depends mainly 
on the monovalent ion and on the temperature. Each trivalent ion is 
surrounded by a nearly regular octahedron of six water molecules. This 
octahedron has its cubic axes slightly rotated with respect to the crystal 
cubic axes and a small trigonal distortion is present along a crystal body 
diagonal. The four different body diagonals provide the four ions in the 
unit cell. 

The 7 orbital levels of Cr+++ in the lowest multiplet are split up in a 
singlet level (J%) which lies lowest, a triplet (J'5) which is & 104 em- 
higher and another triplet (J4) at about 1.8104cem-. For further 
theory of the orbital and the magnetic levels of Cr, reference is made to 
a paper of VAN VLECK [3]. 

The spin Hamiltonian has the form 


(3.02) H 5 = PHS + D[S2—1/38-(S +1)]. 


The g value is nearly isotropic and equal to 1.98. The value of D depends 
on temperature and dilution. BLEANEY [4] observed by paramagnetic 
resonance a gradual transition in the crystal structure below 160° K. 
Two complexes are formed, one with D;=0.135 em- and the second with 
Dy=90.075 em at 20°K. The value of b/C equals 0.75 x 106 92 for the 
concentrated salt. 


[3.2.2] Undiluted CrK alum; H parallel (111) direction. 
The first experiments were carried out with the magnetic field parallel to 
the (111) direction at a frequency of 9400 MHz, using the reflection 
method. For this direction of H one of the four ions in the unit cell has 
its symmetry axis parallel to H, while for the other ions it makes an angle 
of about 70° with H. The resonance spectrum is therefore very much 
spread out. For the concentrated salt the main resonances occur at fields 
of about 3340 9, 4760 6 and 1990 o. 
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The largest resonance at 3340 9 was first saturated and it was found 
that the saturation graphs were curved instead of straight. However, 
at very high and very low powers these graphs were practically straight 
lines. Some examples can be found in fig. 3.1, where only a few of the 
measured points could be plotted due to typographical difficulties. In 
table I the averaged coordinates of the graphs are given, converted to 
the values as would have been found for one gram ion of Crt+++, together 
with the magnetically absorbed powers at minimum and maximum Loe 
The unit of Peis Wate .and of lP,, and n Watt/mole. 

The values of P.w|Pm at P.4=9 are extrapolated and therefore some- 
what uncertain. Their relative values differ in some cases from the expected 
proportionality with temperature. This is due partly to the observed 
experimental conditions. The values of the maximum P,,, are given 
between brackets below the corresponding P,,,/Pm values. Doubtful values 
of 7 are put in brackets. 


TABLE I 
CrK alum H // (111) 


Coord. 
105 x Pext/Pm n= 205 
este, 105 
Pest min max (Pe Ene 
0 0.020 0.050 0.100 0.200 0.400 max Ww 
1.39 36 68 112 148 198 339 WA (oh7f— AUIS GSE QSL tas} 
1.49 36 68 102 158 246 376 626, LO 73) 12899159 0217 225 
(0.800) 
1.78 62 94 1384 192 236 400 10 68 108 147 0.10 192 
1.94 50 72 102 142 16 78 0.09 134 
2.52 158 194 248 300 382 472 586 10 67 170 423 0.12 315 
(1.000) 
2.55 90 122 156 196 249 300 2.5 65 154 480 0.08 180 
3.00 122 150 154 182 228 292 416 18 68 385 515 0.28 250 
(1.600) 
3.03 92 118 150 188 248 14 (75) 81 0.05 150 
3.15 130 146 168 200 274 336 15 68 165 340 0.25 305 
3.22 72 96 132 184 268 358 23 718 25) salon O18 260 
33620 14618259220) 260° 316 400 14 A775 716 0.12) 2380 
(1.900) 
3.16.) U6) 136 7172) 196) 246 7 300 372 5 68 350 920 0.16 2265 
(1.300) 
4.20 124 144 170 208 266 342 400 14 67 £200 (785) 0.15 237 
(0.800) 
1.91 257 
2.28 300- 
380 
2.40 380 


2.71 490 
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[e) aa 400 800 41200 1600mW 
Pext 
Fig. 3.1. Saturation graphs of CrK alum; H // (111) 
A T=2.55°K O T=1,49°K 


a — 3,365 Is 
In fig. 3.2 the values of 7s and 7, are plotted as a function of tempera- 


ture, on a double logarithmic scale and give the results: 
4s=68 W/mole, and yy=77xT!-9 W/mole for 1.39° K<7'<4.2° K 
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At 20.4°K the average over several measurements gave a value of 
7 = 530 W/mole, which is lower than 7; at 4° K, so that apparently very 
high powers have to be applied to reach the final straight part of the 
saturation graphs. 

For the case of internal spin equilibrium the value of the spin lattice 
relaxation time becomes at H =3340 0 according to (1.89) 


2 


Ty = re u ms.deg. 


The saturation of two other resonance lines at 2.55° K resulted in 


7=20 W/mole for H=1990 » 
and 
7=29 W/mole for H=4760 g. 


These values were difficult to determine and are rather doubtful. 

The experiments on the main resonance line were repeated by using 
the transmission method. Essentially the same behaviour of the saturation 
graphs and of 7s and ny was found. 

The values of 7,,, for H = 3340 @ as derived from the relaxation measure- 
ments on powdered specimens of CrK alum by Kramers, Bist and 
GoRTER [5] have also been plotted in fig. 3.2. They found a temperature 
dependence which is practically the same as of our ny while the absolute 
values are also very close. The dependence of t on the magnitude of H 
as found by them is according to the Brons—Van VuEcK formula (1.91) 
with p=0.17 and 0.24 at 3.00° K and 2.51° K respectivily. This holds 
only at lower fields, as at higher ones 7 starts to decrease. 

For not too pure a sample a distribution of relaxation times had to 
be accepted to explain the observed disagreement of the frequency de- 
pendence of the susceptibility with the Casimmr—Du Pr& formulas (1.102). 

For this impure sample the values of t,, were about one fifth of those 
of the pure ones, but the field and temperature dependence was practically 
the same. 

At the higher temperatures measurements have been made by BROER 
[6], by Disxstra, Gorter and Van PaEMEL [7] and by De Vriver [8]. 
The dispersion measurements by Brozr at liquid air temperatures have 
a spreading from sample to sample of about 30 % and when H is varied, 
11 follows the BRons-VAan VuEcK formula with p=0.25. The values of 
» for the purest sample are: 


177° K= 7.7 X 104 W/mole 
n90° K= 1.2 X 105 W/mole. 


The phase transition in the crystal structure near these temperatures will 
influence 7. 
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De Vrijer’s absorption measurements gave: 
n14,3° K= 3.2 X 108 W/mole 
n20.4° K= 3.3 X 103 W/mole 


but due to an uncertainty in the corrections these values are not too 
reliable (possibly a factor 2 higher). To give a better survey these measure- 
ments have been plotted in fig. 3.3. One has to bear in mind that our 
measurements are performed on single crystals, while those of the other 
workers relate to powders. 
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Fig. 3.3. 7 as a function of 7 for CrK alum 
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[8.2.3] Undiluted CrK alum; H // (100) direction. For this 
direction of the magnetic field all ions are magnetically equivalent and 
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their symmetry axis makes an angle of about 55° with H. The spectrum 
consists mainly of three peaks at about 3300 0, 1630 6 and 990 a. At all 
three resonances saturation measurements were performed by the reflec- 
tion as well as by the transmission method. Initially the saturation graphs 
were thought to be straight lines, but in later experiments extending the 
range of P.,, to the lower side all three lines proved to be curved. The 
coordinates of some of them with the values of 7 are given in table IT. 
nf is plotted as a function of temperature in fig. 3.4. 
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Fig. 3.4. 7 as a function of 7 for CrK alum; H // (100) 
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To obtain reasonable values of ns, still lower powers ought to be applied. 
At the minimal power applied their values were of the order of 16 W/mole 
for H=3300 9 up to 4° K. In the liquid helium range the temperature 
dependence of yy is about 7-88 for the two higher field resonances. 
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TABLE II 
KCr(SO4)2.12H20; +H // (100) 


i 


fe on 108+ Pext/Pm 

gaa 7 Pr tee coord. 

P.. i iiaieoainel 7 ; + min ™ max " 1 = Qn 
ext 


1.38 174 234 270 355 477 666 1005 1093 0.23 14.2 32.9 45.3 4.5/372 


(36) 
48 
1.54 211 279 321 388 501 688 1020 1079 0.23 9.9 32.4 —  0.9/277 
(35) 5 
1.82 259 317 363 440 564 763 1099 1134 0.23 16.3 30 54 6.6/524 
(34) 
66 
2.25 300 333 363 415 494 616 — 862 14 22.2 37.2 — 4.5/426 
(32) 74 
2.99 410 463 503 559 640 789 — 1042 0.24 24.7 30.7 63 5.1/581 
(32) 
120 


4.25 448 492 509 560 625 734 947 1044 0.11 16.7 35.4 104 1.1/497 


1630 1.38 78 104 125 156 209 — — 245 O11 5.5 45 8.0 - 
(11) 

1630 4.25 365 394 415 445 485 543 625 661 0.03 3.4 6 21.4 3.0/432 
(40) 26.7 

990 4.25 621 721 803 879 960 1052 — 1071 0.05 0.99 1.6 9.9 2.2/816 


(17) 


The power transfer constant can be represented by: 
np=34 XxT988 W/mole for H=3300 o 
ne= 6.07088 W/mole for H=1630 o 
ng= 2.4X7T9.88 W/mole for H= 9909 


These values of nr when plotted against b/C + H2 on a double logarithmic 
scale lie on a straight line: 


ny =1.1(b/C + H2)1-4, 


The power transfer might perhaps be due to two or more parallel transfer 
coefficients, each with a different temperature dependence. We tried a 
compound curve n= 42.5+4.5 71-9 W/mole as for the (111) direction 
the yy has a 71.9 law. This curve has been drawn in fig. 3.4 and seems 
to fit the experimental points of the line at H—=3300 9 rather well. 
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[3.2.4] CrK alum diluted with AIK alum. To find the 
dependence of 7 on magnetic dilution, two samples have been measured 
with Cr : Al ratios of 1 : 2.28 and 1 : 40. The direction of H was parallel 
(iab 1): 

ESCHENFELDER and WEIDNER [9] obtained values for a 1 : 10 diluted 
sample, namely 71;10=66 W/mole independent of temperature between 
2° K and 4° K. They saturated the spin system to a very high degree 
(22 times), but their data were only obtained from selected straight 
saturation graphs, and no information was given about the curvature. 

The results of our measurements can be found in table IIT and in 
fig. 3.5. In this figure we have also drawn 71:19 of ESCHENFELDER and 
WEIDNER. 

The Cr*+*+* concentrations were determined by the polarographic 
method *) [29]. 


TABLE III 
(Cr/Al)K alum; H// (111) 


Cr Ale 528 Cr: Al = 1:40 
min max min max 
1.89 133 Dal 110 2.00 58.2 35 51 
Pies 218 5.0 init! | 2.99 58.9 20 43 
3.47 377 15 141 AD, 55 14 47 
4,22 445 26 162 14.6 w117 Dill 74 
20.0 123 16 4] 


Ae Dilution Or Al=a1-= 2,23: 

All saturation graphs are curved at low input powers. The curved part 
was small, so that reasonable values of the starting slope could not be 
determined. The main part of the graphs is straight, with reciprocal slopes 
as given in table III. 

The absorbed power at which the curves become straight is about 
15 W/mole which is much less than for the undiluted substance with 
H // (111). We find: 

71:2.28= 91 x 71.58 W/mole. 


By the relaxation method Van DER Maret e.a. [10] find practically the 
same temperature dependence (7'1:8) for a sample diluted to the same 
ratio. The absolute value of 7,,, as found by them is anomalous, as they 
are smaller than for any other concentration used at most of the magnetic 


fields and temperatures. 


*) We are very much indebted to Dr W. L. Groeneveld and A. P. Zuur of the 
Chem. Dept. of the Leiden Univ. for the use of their equipment and their advice 


on these matters. 


4 10 2.5 K 


V Cr:Al= 1:2 relaxation measurements 
—— Cr:Al 1:10 ESCHENFELDER and WEIDNER. 


b) Dilution Cr: Al=1 : 40. 


The saturation graphs are straight lines. The molar power transfer 
constant is found to be (ef. fig. 3.5) 


m1:40= 58 W/mole for 2° K <7 <4° K. 


At two hydrogen temperatures measurements were made (fig. 3.5). 
The highest dilution measured by VAN DER Maret is 1 : 20 for which 
he found 7,.,=62 W/mole at 2.9°K and H=3375 9 and a 7'1-6 law at 


low magnetic fields. At higher magnetic fields this dependence becomes 
somewhat smaller. 


[3.2.5] The conclusion drawn from the above mentioned experi- 
ments is that for CrK alum nz equals 7,., whilst 75 behaves as theoretically 
expected. From the measurements with H//(111) and H//(100) an average 
value of ny is derived, very close to the value of 7, for powders. From 
fig. 3.3 it is seen that the values of yy extrapolate nicely from the liquid 
helium range to liquid air temperatures with a 71-9 law. Extrapolation 
of 7s would result in a 7’8-4 dependence. Both powers of 7’ are much lower 
than the theoretically expected value of 6 for a purely indirect process. 
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[3.3] Results on CuK and CuNH, tutton salt 


The relaxation behaviour of CuK tutton salt has been studied as a 
function of temperature, direction of the static magnetic field with respect 
to the crystal axes and preliminarily of the microwave frequency. A 
rigorous test was set to the apparatus by making the measurements 
over a wide range of applied powers. 


[3.3.1] The tutton salts form an isomorphous series with the 
general formula M++M2+(S'VOs)26H20 where M+ and M++ are mono- 
valent and divalent cations resp., and S'V=S or Se. The crystal is mono- 
clinic (Turron) [11]. Each divalent cation is in a crystal electric field 
of nearly tetragonal symmetry caused by the surrounding distorted octahe- 
dron of six water molecules. There are two differently oriented magnetic 
complexes M++(6H20), the ac-plane being a plane of reflection. 

The angle of the magnetic z-axis with the ac-plane is called « and the 
angle between the projection of z on the ac-plane (K, axis) and c-axis 
is called y (cf. fig. 2, ref. [12]). The angle is measured in that sense .so as 
to make the angle between a and c obtuse. 

The magnetic data of the CuK and CuNH, tutton salts are assembled 
in table IV. The g values and « and yw have been found to be temperature 
dependent [13]. This is interpreted as due to a change in the crystal field 
which brings a different orbital level to the bottom [14]. PoLpER’s [15] 
calculations on the g-factor are in fair agreement with experiments. 


TABLE IV 
Units of CT2/R are (degrees)? x 10-4. 


Ge gy ge «& wp COmT?/R COnucT?/R OaipT?/R CoxT?/R 


CuKe(SO4)2 6H20 —«-2..31 2.07 2.25 41° 105° ~—- 6.0 Ld ine 3.6 
CuNH,(SO.)2 6H2O 2.32 2.09 2.25 39° 77° 8.7 eal 3 6.3 


He assumed a tetragonal field of the same order as the cubic one and showed 
that the sequence of the orbital levels starting from the lowest one, must 
be I3, 11, I; and I's and that the energy level separation of [4 and J3 
is 15400 cm-1 and that of I’; and J3 is 26600 cm-!. It has been shown, 
however, that a rhombic term in the crystal field is present [13]. The 
spin specific heat contains a nuclear, dipolar and an exchange part and 
has been analysed by Brenziz and Cooke [16]. 


[3.3.2] CuKe(SO4)2 6H20. Some typical saturation graphs are 
shown in fig. 3.6, and the coordinates of some of the graphs are given in 
table V. The plotted curves contained about three times as many points 
which could not be inserted in the drawing. The reproducibility of the 
apparatus is demonstrated here by the fact that some points of repeated 
series are three fold coinciding. The curvature of the graphs is seen (from 
fig. 3.6 and table V) to depend on the direction of the magnetic field. 
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Fig. 3.6. Saturation graphs of CuK tutton salt. 
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The minimum P,,, used was of the order of 5 to 10 uW, but to obtain 
values of 7; for all temperatures and directions of H, still smaller powers 
ought to be applied. In fig. 3.7 the values of 7 obtained for H parallel to 
the crystal a-axis (magnetic K1) is plotted as a function of the tempera- 
ture. Two sets of measurements were performed, the earlier ones with 
the reflection method using very high powers, while with the later ones, 
by the transmission method, also low powers were used. This may explain 
the difference between the results obtained by the two methods. The 
discrepancy is not too serious, however, as the measurements were made 
on different crystals and as the results of VAN DER MAREL e.a. indicate 
an influence of the crystal condition. His results for 7,,, obtained with 
the relaxation method, for a freshly grown powder and a slightly deteri- 
orated one, are also shown in fig. 3.7. For this last sample 7, is about 
1.6 times as large as for the first one. 

To check the influence of cooling of the crystal by liquid or by gaseous 
helium some runs were made with the liquid kept out of the cavity the 
crystals being glued to a cavity wall. There is a difference of about 20 % 
in the 7 which is of the order of the experimental error and which may 
be due to the use of different crystals. 


TABLE V 
CuK tutton salt 


108- Pext/Pm 


H // (100) > 6 = 90° 


IPs, IP 

Pent i min ty max 

0.1 0.2 0.4 0.8 1.6 ah max 
1.46 318 425 522 688 974 1478 2450 2694 0.93 0.11 413 
(3.6) 
9.08 515 582 650 778 1019 1465 2256 2460 1.49 0.11 1.88 1.4 
(3.5) 
2.51 728 786 844 960 1185 1607 2329 2490 1.81 0.07 2.04 1.4 
(3.5) 
3.48 535 615 685 821 1061 1487 2226 2523 2.5 0.10 262) Wee 
(3.6) — 
2.25 
4.25 614 682 748 881 1134 1542 2295 2486 2D 0.12 Dy LA! 
(3.5) 

Pext Pm Pin 
mW 0 1 2 4 8 16 32 max Ns es "F oes 
oe | a 
92° 794 1079 1337 1819 2636 4069 — 5079 Oars LO sea h.0 3.0 

(2.2) 
90° 556 1635 2323 3351 4920 7549 1.82104 1.76104 0.39 eA! Ome Ondo 
(4.5) 
112° 603 2047 2959 — 4387 0.13 9.2x10-3 1.05 8 
(3.5) 
140° 524 1270 2016 3224 — — -- 4653 1.05 1.9x10-3 1.74 1.4 


(6.5) 
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Fig. 3.7. 7 as a function of T for CuK tutton salt; H // Ky 


nf [st series gas cooled 

nr [st series liquid cooled 
nr 2nd series liquid cooled 
Nrer. powder 

Nre: deteriorated powder. 
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The temperature dependence of x in the liquid helium range as found 
by us, with H // (100) direction, is 7°48 which is the same as for 7,,. This 
temperature dependence of 7 seems to change to 71-9 above 4° K. 

The dependence of yz on the direction of the magnetic field with respect 
to the crystal axis, was measured in the following way (transmission 
method). With a fixed applied power the values of x2 and 6 were measured 
as a function of the angle @ between H and the crystal axes, H being 
rotated in a plane perpendicular to the magnetic r.f. field. This was 
performed at a very high value of P,,, and at an intermediate one. This 
last level was taken so as to be well in the final straight part of the 
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saturation graphs. From these two measurements nr aS a function of 6 
can be obtained. This procedure was repeated several times. As a check 
a complete saturation graph was made for some angles. 

For H rotated in the ab-plane (K,K3) the effective g-factor (from the 
position of the resonance line due to both ions) is seen to have the required 
symmetry (fig. 3.8). Due to the facts that the line widths for the two 
inequivalent ions are not the same and that 2« 490°, the effective g -factor 
is not constant. 7¢ shows about the same symmetry as ¥,-, the minimum, 
however, appears to be slightly displaced (about 15°) with respect to the 
tetragonal axis. The variation in np due to both ions is of the order of 
15 %. However for each ion separately this will be larger. 

For H rotating in the ac-plane (K1K2), ns changes by a factor 4 (cf. 
fig. 3.9), but the angles at which the extremes occur do not correspond 
with the magnetic axis, the displacement being about 30 to 40 degrees. The 
direction of maximum 7; corresponds to the absolute minimum of the line 
width, and the direction of the cutting plane of the crystal (see also [13)). 

The angular dependence of nF is seen to be in disaccord with Nagaoka’s 
calculations (cf. [1.5], [17]) as he expected 7 to become very small in the 
direction of the magnetic axes. The few data on the angular dependence 
of 7s seem to suggest that this one has the same variation with 6 as the 
g-value curve in this plane. 

Preliminary measurements on the frequency dependence of 7 indicated 
this to be an increasing function when changing w/2x from 8300 MHz to 
10 500 MHz. 

For this salt it may be concluded that ny and 7, are equal. 


[3.3.4] Cu(NHa4)2(SO4)2 6H2O. To examine the influence of the 
exchange interaction on the relaxation behaviour, experiments were per- 
formed on CuNH, tutton salt, where this interaction is nearly twice as high 
as in CuK tutton salt, the dipole dipole interaction being equal (cf. table IV). 
The saturation graphs are curved and shown in fig. 3.10, where only a 
few of the measured points could be inserted. 

The magnetic field was parallel to the crystal a-axis. At 1.46° K and 
1.83° K the cavity Qo was differing from the value for the other tempera- 
tures, reasons why the ordinates of their saturation graphs are relatively 
higher with respect to the other curves, for which the ordinates are 
proportional to 7 for a low value of P,,,. 

The values derived for 7 would have been more satisfactory if the 
range of P,,, had been larger especially on the lower side. The results are 
plotted in fig. 3.11. Little value has to be attached to the results of ns. 

At temperatures below 2.4° K ny is seen to have the value of 2.26 W/mole, 
independent of temperature. This value is about a factor two lower than 
for CuK tutton salt at 1.4° K, whilst 7, is less dependent on temperature. 
The different orientation of the magnetic axes makes comparison difficult, 


however. 
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Fig. 3.10. Saturation graphs for CuNH, tutton salt; H // (100). 
Temperatures next to graphs. 


GIORDMAINE e.a. [18] deduced, from decay and steady state saturation 
measurements on a set of crystals diluted with the Zn salt (concentrations 
Cu : Zn up to 1 : 10), relaxation times of 20 s. from decay and 2 s. from 
steady state measurements at 1.3° K and near 9000 MHz. The steady 
state saturation value has to be multiplied by 8 when homogeneous 
saturation takes place as was indeed observed by him. The corresponding 
7 would be about 0.02 W/mole, increasing with concentration. This value 
is two orders of magnitude smaller than our 77. 


(To be continued) 
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ON THE POWER TRANSFER BETWEEN PARAMAGNETIC SPINS 
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(Communicated by Prof. C. J. Gorter at the meeting of June 27, 1959) 


[3.4] Results on CuSO, 5H20 


The exchange interaction in this salt is larger than the Zeeman 
energy used, and one therefore expects the relaxation times to be very 
short and of the order of t2 (see [1.4.3]). 
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Fig. 3.11.» as a function of T for CuNH, tutton salt; H // (100) 
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[3.4.1] The crystal structure has been studied by Brnvers and 
Lipson [19] and is found to be triclinic. The cupric ion is in an electric 
field of nearly tetragonal symmetry, caused by four HO molecules in a 
square at a distance of 2 A and two oxygen atoms belonging to the SO, 
groups further away. The unit cell contains two Cut++ ions having an angle 
of 82 degrees between their tetragonal axes, but resonance and suscepti- 
bility measurements indicate, that this angle is closer to 90°, making the 
crystal magnetically uniaxial. The cupric ion is in a D state, which is 
split up in a cubic field in a triplet and a doublet, the latter one lying 
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lowest. In the tetragonal field the doublet is split up with an estimated 
separation of the order of 12000 em-! [15]. 

The mininum susceptibility axis LZ; makes angles of 154°, 64°, 51° 
with the triclinic axes (a,b,c). The resonance spectrum has been in- 
vestigated by BaGGULEY and GriFFITHs [20] at various frequencies, they 
found gz=2.46 and gz, y=2.08. 


TABLE VI 
rrr 
sol. 1 sol. 2. interaction with 2 ions of 
Ji/k 02935 0.97° the same kind at 5.97 A. 
Je/k —0.36° —0.23° the same kind at 6.12 A. 
J 3/k 0.085° —0.085° a different kind at 5.19 A. 


In an analysis Pryce [21] suggests that one has to distinguish the 
neighbours into three pairs, each pair having a different value for the 
exchange constant J. The possible solutions for J are given in table VI in 
units of °K. They are obtained by considering the line width, ranging 
from 50 to 220 @, the splitting of the line at higher Larmor frequencies in 
certain directions of H, the specific heat Cy =0.75 R/T? (degr.)? mole [22] 
and the Curie temperature of 6=—0.65 °K. AsHMEAD [23] has found 
anomalies in the specific heat near 1°K and 0.01 °K. 


[3.4.2] The saturation graphs were all curved, but the range in 
which this curvature occurred, could not be reproduced very well for 
different helium runs and for different crystals. However, the values of 
ns reproduced rather well. They are plotted in fig. 3.12 together with some 
values of 7 and are seen to decrease faster than 1/7’. The direction of H 
was close to the projection of Z; on the (110)-plane (5° off). We have 
also plotted 1/7; vs T on a linear scale in fig. 3.13. The points lie very 
close to 


1/ns= 1.4 x 10-2(T' — 0.94) (W/mole)-1. 


This resembles a kind of Curie-Weiss law with a value of 0= + 0.94° K, 
which is close to the maximum of the specific heat anomaly at 1° K and 
the temperature at which one might expect Ji to freeze out. Just above 
such a transition large fluctuations due to short range ordering occur 
in the local fields which might cause the increase of 7 near this transition. 
To check whether the relation 7” « yo really holds at low values of P.,, 
(cf. [1.7.6]), a plot was made of P,,,/Pm versus T at P...=5 mW the 
lowest P,,, common to all curves considered. The points lie close to the 
straight line A(7'+0.5). The agreement is satisfactory. 

For the best series we plotted P,,,/Pm against /P.., as a trial. The 
points are very close to a straight line and resemble the curves of 


CuK2Cl, 2H20 to be presented later in fig. 3.14. 


Fig. 3.12. 7 as a function of T for CuSO4 5H2O 
eT, O 


ae 


Fig. 3.13. CuSO4 5H20; I/n. vs. 7. 
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Fig. 3.14. Saturation graphs of CuK2Cl, 2H20. 
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[3.5] Results on CuKeCl, 2H20 and Cu(NH4)2Cly 2H20 


Measurements on these salts have been carried out partly because the 
exchange interaction is intermediate between that of copper sulfate and 
of CuK tutton salt, partly as preliminary measurements for a double 
resonance experiment to be discussed later. 


[3.5.1] The crystal structure has been studied by CuroBak [24] 
and Wycxorr [25]. The crystal symmetry of these two salts is tetragonal. 
A unit cell contains two cupric ions each in a tetragonal field due to four 
Cl- ions in a rhomb lying in the crystal ab-plane and two HzO molecules 
along the normal (c-axis). For the two Cut++ ions the rhomb is rotated 
(around c-axis), so as to make the ac-plane a plane of reflection. For 
CuKeCl,2H20 the g values are gy,-=2.06, gz=2.105 and g,=2.30. The 
anisotropic line width AH ranges from 110 to 250 o. In the (001) direction 
AH=1159. The average exchange constant J is estimated to be 
0.0931 em-! = 2792 MHz [26], but the values of the various exchange 
interactions present have not yet been analysed. For Cu(NH4)sCly 
2H2O these values are somewhat different. 


[3.5.2] CuKeClya 2H20. Atsome temperatures the saturation graphs 
have been measured over a large range of P,,,, their coordinates being 
given in table VII. The direction of H was at 50° to the c-axis in the 
plane bisecting the two a-axes. Some of the measured graphs have been 
drawn in fig. 3.14. In this figure are also plotted the values of P,,,/Pm 
as a function of /P..4. These last points are very close to a straight line 
in the central part of the saturation graphs which are therefore close to 
parabolas with their start and final parts straightened. This was the case for 
all graphs of table VII. The slope of the straight line is given in the last 
column, together with the ordinate at P.,,=0. The values of P.,,/Pm as 
a function of 7 at P..,=3 mW lie on a straight line through the origin. 

To determine more fully 7s and »y more experiments were performed 


at only low or high values of P,,,. For different crystals the curved part 


TABLE VIL 


CukKsCl, 2H2O 


108. Pext/Pm 


1.87 511 552 586 640 716 830 982 — 1078 3.4 50.5 42 244 
(45) 

3.22 667 737 773 818 888 989 1137 1310 1474 2.6 33 73 220 
(100) 

4.25 — — 1160 1191 1245 1323 1399 1476 1584 1.7 23.2 71 660 


(110) 


(ihe 


was not completely in the same range of P.., like for CuSO, 5H20, but 
ns and ny were reproducible. They are plotted in fig. 3.15 and fig. 3.16. 

In the helium range we find when plotted on a double logarithmic 
scale ns ~ 95/T W/mole. But, when making a linear plot of ys vs 7, a 
slightly better fit is obtained by the formula: 


IZ 


Uh = TOA W/mole. 


The values of 77 suggest a discontinuity at the 2 point of helium, (although. 
we do not know if this is a reality) with a 73 or 74 dependence above 
this point. 

For H // c-axis the value of 7 was determined at a few temperatures, 
and found to be four times smaller than for the foregoing direction in- 
dicating a strong anisotropy of 7. VAN DEN BROEK e.a. [27] measured 
Ne On a large single crystal. Their values for H=3300 6 and parallel to 
the c-axis have been inserted in fig. 3.15 showing a dependence and the 
discontinuity at the A point. So for this salt we find the relation 7,,. ~ nf 
to be fulfilled when taking the anisotropy into account. 
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Fig. 3.16. 7. as a function of T for CuKe Cl42H20 


TABLE VIII 
Mn(NH4)2(SO4)2* 6H20 


108+ Pext/Pm 


Pi Ps coord. 
Pext min Be max . n = 2ns 
mW 0 O42 O88, 915 3.0 6 12 max 
T°K | 
37.6 
1.33 99.8 149 182 229 306 427 617 914 0.88 8.7 25.2 — _ 1.0/196 
(23) 34.4 
14.7 
159 126 162 184 214 268 360 528 786 0.23 — 245 36 0,5/167 
(21) 13.6 
2.06 181 196 208 216 260 $21 441 461 O21 21.7 27.9 50.8 - 
(13) 
2.65 228 253 260 269 286 318 — 368 0.42 81 29.4 91 0.2/243 
(10.5) 
3.31 283 300 311 324 342 371 — 417 0.39 11.9 29.4 99 0.75/309 
(10.5) 
4.25 


333 361 370 379 392 415 458 495 0.09 10.6 38.6 136.3 0.3/357 
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[3.5.3] Preliminary results on Cu(NH4)2Cl4 2H2O indicate ny 
to be about the same as for the potassium salt. Also the temperature 
dependence has the same character. At 17° K it was found that 7 « 400 
W/mole which seems rather low when extrapolating n,. 


[3.6] Results on Mn(NH4)2(SO.)2 6H2O . 


The manganese ion has a spin of 5/2 and is an S state ion. One has 
to go to a high order in perturbation theory to find the interaction between 
spins and lattice, and therefore expects to find small values of 7. 

One of the reasons for investigating this salt was the near independence 
of the temperature behaviour of 7 on magnetic dilution as found by 
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Fig. 3.17. 7 as a function of 7 for MnNH, tutton salt 
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VAN DER Maren e.a. [10]. Only the concentrated salt has yet been 
investigated, but diluted crystals will be measured. The crystal structure 
of the tutton salts has been discussed in [3.3.1]. The angles mentioned 
there are ~«—32° and y=58°. Due to the fact that S=5/2 one has to 
add terms of the fourth power in the spin operators to the spin Hamil- 
tonian. The g-value is 2.00. 

The saturation graphs are curved and the coordinates are given in 
table VIII. The resulting values of 7 are plotted in fig. 3.17 together 
with the values found by VAN DER Maret. His sample contained imputri- 
ties, probably of a physical nature, which increased the value of 7,,). 
The value he finds for b/C is also 16 % smaller than the ones of BROER [6] 
and of Bist [28]. The results of the relaxation measurements for H = 33759 
are also inserted in fig. 3.17. VAN DER MAREL e.a.’s values at hydrogen 
temperatures are seen not to correspond with the extrapolation of their 
helium measurements. Also the physical purity has been found to affect 
strongly the shape of the 7(@) vs w curves, which do not always correspond 
with true Casimir—Du Pré curves. The measurements of B1sL corresponded 
satisfactorily with Casimir—Du Pré curves and 7 follows the Brons—Van 
Vleck formula (1.90) at low values of H with p=0.43. We find: 


np=19.5 T1-4 W/mole 


which also coincides with van der Marel’s values at hydrogen temperatures. 
The values of 7; are somewhat uncertain as smaller P,,,’s are necessary 
for its exact determination. 


[3.7] Results on Co(NH4)2(SO4)2 6H2O 


The strange results on this salt reported in an article by VAN DEN BROEK, 
VAN DER Maret and GorTER [12] tempted us to investigate this salt, 
which is especially different from the other salts in its dependence of t 
on the magnetic field. t is found to decrease initially with magnetic field 
and to have, for the concentrated salt, a 7-1-0 temperature dependence, 
while at higher magnetic fields it increases with H and has a 7-2-4 
dependence. For the details we refer to the above-mentioned paper. 

The crystal structure of the tutton salts has been discussed in 
[3.3.1]. The Cot+ ion has in a cubie field an orbital triplet lowest. This 
triplet with its fourfold spin degeneracy is split into Kramer’s doublets 
by fields of lower symmetry. The effective spin S’ equals 1/2. 

The data are y= +137°, «=34°, and the g-values along the principle 
susceptibility axes are: g1=5.70, go=3.06, g3 = 4.36. 

The measurements were made with H // Ky axis. The saturation graphs 
are curved as is evident from fig. 3.18, the coordinates are given in table IX. 
Only a few of the measured points are inserted in the drawing due to 
typographical difficulties. The extrapolated ordinates at P.4=0 are 
except for the graph at 4.2°K proportional to 7’. The peau values 
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TABLE IX 


Co(NHa4)2(SOa)2 ° 6H2O 


108+ Pext/Pm 


3.30 


4.25 


316 384 427 


394 454 488 


551 604 639 


590 628 655 


724 755 780 


757 810 846 


200 


571 671 858 1010 0.57 21.8 
(99) 
620 734 907 1106 0.50 22.8 
(110) 
760 870 1040 1256 0.71 34.5 
(108) 
749 832 968 1139 0.17 53.6 
(115) 
904 1022 1232 1331 1.36 54 
(80) 
960 1074 1256 1314 0.39 85 
(75) 
Oo 
= Oo 
ie. = _| 


2 4 6 Ok 


161 


182 


Fig. 3.19. 7 as a function of 7 for CoNH, tutton salt; H // Ky 
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of ms and ny are plotted as a function of temperature in fig. 3.19 together 
with the values found on powdered samples by VaN DEN BRoEK e.a. 
Their values correspond very closely to our 7s. This is the first example 
where the value 7s, derived from the starting slope, equals 77... The values 
of ny decrease slightly with temperature. 

As the coupling of the ion with the lattice is strongly anisotropic it 
would be interesting to determine the anisotropy of 7. 


(To be continued) 


PHYSICS 


ON THE POWER TRANSFER BETWEEN PARAMAGNETIC SPINS 
AND CRYSTAL LATTICE. IIIc 


BY 


B. BOLGER 


(Communicated by Prof. C. J. Gorter at the meeting of June 27, 1959) 


B. Results obtained by pulse measurements and relaxation rates determined 
by experiments on masers 


[3.8] Results on KgCr(CN)g/K3Co(CN)¢ 


Measurements were performed on this salt for different chromium 
concentrations, among which sample A containing 0.035 % Cr to Co and 
sample B with 0.1 % Cr to Co. The measurements described here were 
performed with a signal frequency of 1420 MHz and a pump frequency 
of 3850 MHz. 


[3.8.1] Crystals of KsCr(CN)g¢ diluted by the isomorphous Co 
salt have been studied by Baker, BLEANEY and Bowers [30] and 
recently by many others because of their application in masers. The 
zero field splitting of the Cr+++ion with S =3/2 corresponds with a frequency 
of 5700 MHz. There are two chromium ions in a unit cell, with their 
magnetic axes rotated about the c-axis by + 6° with respect to the 
crystal axes. In fig. 1.1 the energy levels were shown as a function of 
the magnetic field. In fig. 1.2 the loci are given of some of the absorption 
lines in a plane of the H space, measured at frequencies of 1420, 3850 
and 8600 MHz; the magnetic vector H is rotated in the ac-plane of the 
crystal. The energy levels giving rise to the transitions are shown beside 
each locus, the levels being numbered in order of increasing energy. If the 
1-3 transition is saturated at 3850 MHz there is seen to be a possibility of 
maser action at 1420 MHz with the 1-2 transition when the two resonances 
coincide, which happens for two field conditions in each quadrant, and 
with the 2-3 transition at another two field conditions. The operating points 
will be denoted I, I, HI and IV in order of increasing magnetic field. 
It is also possible to saturate the 2-4 transitions at 8600 MHz and use 
the 2-3 transition for the signal frequency of 1420 MHz. 

The line width of the 1-2 transition at 1420 MHz for both samples 
was of the order of 8 9. Especially for sample A this width is much larger 
than one expects on the basis of homogeneous line broadening (& 0.8 9) 
(cf. [1.4.4]). One may therefore conclude that the lines for sample A are 
inhomogeneously broadened. 
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The matrix elements of the different magnetic dipole transitions as a 
function of the orientation and magnitude of H have been computed 
by HowarrtH [31]. 

The chemical analysis of the chromium content was performed by 
the method described by Dan [32], and by rontgen fluorescence spectro- 
graphy *. 


[3.8.2] The influence of a modulated pump power on the 
signal. One has to work with pump and signal frequencies, which are 
fixed by the cavity resonances. Finding the exact magnitude and direction 
of H for which the pump influences the signal most, e.g. the operating 
points, may therefore present some difficulties. We found a quick method 
to trace the operating points by pulse modulating the pump power at a 
repetition rate lower than the relaxation rates. Near an operating point 
the signal power will also become modulated and the largest modulation 
will occur at its centre. 

With the arrangement of fig. 2.6, using an unmodulated 1420 MHz 
signal, the time behaviour of the signal absorption is displayed on the 
oscillograph. This one, a tektronix type 532, had its time base synchronized 
to the pump modulation, and pictures like those in fig. 3.20 may be 
observed. Photographs were taken from the oscillograph pictures and 
afterwards projected on a screen. Two methods can then be followed to 
determine the time constants t. The first is to measure the deviation of 
the signal from its equilibrium value as a function of time and plot this 
on a single logarithmic scale as shown in fig. 3.21. The second, a much 
faster but less accurate method, is by trying to identify the projected 
pictures with one out of a set of exponential curves with different time 
constants, drawn on a glass plate. One should only know the time scale 
of the projection which is obtained from the projected oscillograph screen 
and its calibrated scale. Both methods to determine t have been used. 
The values of the magnetic field and its direction for the four operating 
points are given in table X. 


on fott 
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Fig. 3.20. Signal absorption as a function of time for a modulated pump power ee 
a) KgCr/Co(CN)e b) Ruby operating point Ry 


*) Tam very much indebted to Drs J. de Vries of the N.Y. Philips’ Gloeilampen- 
fabrieken for performing this analysis. 
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Fig. 3.21. KgCr/Co(CN)s; Time behaviour of x” after a pulse of pump power 


TABLE X 


Data of operating points in KsCr(CN)6/KsCo(CN)« 
Pump 3850 MHz (1-3 tr.); Signal 1420 MHz 


Operating point I Il Iil EV; 
Dignal. transition. ae ee eee as 1-2 2-3 2-3 1-2 
Magnetic field (9) .. . 5 Bice’ eared 480 680 1170 1280 
Z.H to a-axis (in ac- bere Sarr: 11° 25° 14° 3° 
ace fo) tim {U2 ) Pyeng 9 
A Or/Co = 0,088 %.....5.. s = 18 —105  —D.00 —1.3) 
Jeg MOrdk Clo) =A Me be nS 0 0 0 0 
OM Cr/ Co = /Ol0b5 9, 2s et ee 0 0 w—l 
For sample A 
CealUigs 8 ns i, sae Uae vee 4.6 1/3.1 1/2.3 5.7] 
Top (mos.) at 1.4 °K see ink Soe Se ee BAT 26.1 33.4 29.0 


In the last row the values of the rise times of the signal transitions 
Tor Curing the off time of the pump power are presented for sample A 
at 1.4° K. The relaxation time as a function of temperature at operating 
point I was determined for both samples. The results are given in fig. 3.22. 
For both samples we may take rt « 7'-3.3 between 6° K and 20° K_ but 
below this range the temperature dependences differ. For sample A we 
found at liquid helium temperatures 


(3.03) ta=62/T ms. 


For sample B the values of tg in this temperature range can very well be 
represented by 


(3.04) te 2 (gy 


This curve has been drawn in fig. 3.22. 

It was checked whether the pulse length of the pump modulation was 
of any influence on the characteristic times. For sample B at nearly all 
temperatures investigated the pump pulse was varied between 1 ms and 


Fig. 3.22. tj, as a function of 7 for KgCr/Co(CN). 
OO Tog sample A O Tog sample B 


100 ms. The rise time of the signal t,, was found not to depend on the 
pulse length although there is a slight indication of a decrease below 
~~ 4ms pulse length. 

We will now discuss the unmodulated maser experiments in so far as 
they give information about the relaxation properties. 


[3.8.3] The influence of an unmodulated pump power. 
Information about the relaxation rates can be obtained by determining 
the ratio of the signal absorption with P,,,,,=0 to that when P,,.,, is 
very large and by using formula (2.31). For the four operating points 
possible when wp/2x = 3850 MHz and w;/27=1420 MHz these ratios 
are given in table X together with the calculated values of U3/ U2, the 
signal transition being indicated above each column. At some field con- 
ditions for samples B and C zero absorption was obtained for P,.., 
large (cf. table X). Evidently no maser action is possible at those field 
conditions for the Cr concentration used, because of strong cross relaxations. 
From the gain saturation of the maser when large signal powers are 
applied and from the maximum power emitted by the sample more 
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information can be obtained about the relaxation rates. The maximum 
power emitted Piz... provides reliable information about Uos3 as it only 
contains the sample weight but not matrix elements or filling factors 
(see (2.32)). For the operating point no. I of sample A, the only point 
for which amplification was obtained, the measured value was 


Pi2ax = 210-7 W. 
(Qo = 4000, Qe=1000, Qm=—900). 


The calculated value (with (2.32)) for a crystal of 6.45 g, and a Cr/Co 
percentage of 0,035 % as used by us, becomes 


Pie, ,ax = 1-6 X 10-8 X U2g W 


using the ratio of U23/Ui2 given in table X. Comparing the measured and 
calculated values of Po3,,,, results in: 


Uo3= 2.6 s-l and Ujy2=0.6 g-l, 


For the other operating points —wm was for sample A not large enough 
to obtain amplification. For sample B the absorption approached 
asymptotically wm=0 when increasing the pump power. 

A value of Uy.+ U3 can be derived from the saturation of the maser 
gain. As large samples are used in masers, one has to take into account 
the variation of the magnetic h.f. field over the sample size. The resulting 
value will not be too reliable. 


[3.8.4] In-homogeneous saturation at the pump frequency. 
A crucial experiment was performed to test whether a homogeneous 
saturation by the pump is effected even when there are nodes of the 
magnetic field at w, inside the sample. By using the hot phonon theory, 
to be discussed later, Morris, Kyui and StranpBeRG [33] suggested 
that, because of the large mean free paths at low temperatures, phonons 
in a narrow band around w, will transport energy from saturated towards 
unsaturated regions, thus providing a homogeneous spin temperature at 
the pump frequency throughout the crystal. They used this picture to 
explain the large values of —um obtained by them in a maser using a 
ruby and w)/27=24 KMHz. 

Because of the importance of this question in the design of a maser 
and the understanding of the relaxation mechanisms we have carried out 
a simpler experiment than that performed by Morris e.a., where it was 
not possible to measure the emission for a saturating field uniform 
throughout the crystal. 

In figure 3.23 shows a section through a rectangular cavity (cavity no. 3) 
resonant in the TEj93 mode at 3850 MHz. In the cavity a resonant strip 
is placed approximately 4/2 long at 1420 MHz, and a large crystal of 
K3Cr/(CN)¢/KsCo(CN)g with its centre at a point where the magnetic 
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field at the TEj93 mode is a maximum. As indicated on the figure the 
crystal extends about 2/3 of the way towards the magnetic field nodes. 
Under these conditions most of the crystal could be saturated with a 
power of 100 mW in the cavity at 3850 MHz, and the corresponding 
emission was measured (see fig. 3.24). For this cavity the TE\94 resonance 
is at 4100 MHz, with the magnetic field distribution indicated at the 
bottom of figure 3.23. The magnetic field maximum is now near the edge 
of the crystal, and a field node lies within the crystal. Under these con- 


K,Co (CN), Crystal | | 


h=o 


Fig. 3.23. Cavity for inhomogeneous pumping experiment 


A TE103 mode 3850 MHz 
B TEi04 mode 4100 MHz 
C coupling slot 1420 MHz 
D coupling loop 3850 MHz 
S strip 4/2 at 1420 MHz 
h = 0 nodes of pumping field 
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P, 
Fig. 3.24. Signal absorption as a function of pump power, Py 


[J pumped by 7'H103 mode (no pump node in crystal) 
© pumped by T'H104 mode (pump node in crystal) 
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ditions there was no net emission with a pump power of 100 mW in the 
cavity, at which level the wm was close to an asymptotic value, but still 
decreasing with increasing P,,,, a8 can be seen in fig. 3.24. 

An auxiliary experiment with another cavity confirmed that the 
emission from a small, uniformly saturated crystal was comparable for 
saturating frequencies of 3850 and 4100 MHz. These experiments suggest 
that, at least for KgCr(CN)¢/K3Co(CN). and with these pump frequencies, 
it is not possible to saturate a crystal throughout its volume with a non- 
uniform field distribution. The experiment by Morris e.a. might be 
explained by realizing that for large pump powers the regions where the 
pumping is insufficient become relatively small, which is demonstrated 
by the behaviour of wmioa (cf. fig. 3.24). 


[3.9] Results on synthetic ruby 


Synthetic ruby, AlgOs with a small amount.Cr203, has proved to be a 
suitable maser material. 

The paramagnetic resonance spectrum of the Cr+++ ion in this sub- 
stance has been studied by various investigators [34] [35] [386]. The zero 
field splitting corresponds with a frequency of 11593 MHz with no rhombic 
component being present and with an isotropic g-value of 2.00. In 
fig. 3.25 some of the loci of the transitions in H space are shown (H in 
ac-plane). 

For the 1-2 transition at 1430 MHz it is possible to use the 1-3 transition 
at 10.7 KMHz as a pump frequency. An operating point occurs for 
H=430 9 and at 23° to the c-axis of the crystal, which point will be 
denoted by R,. Another operating point investigated, Ry, was obtained 


c-axis 


Fig. 3.25. Rotation diagram of ruby, symmetry axis is c-axis. Corresponding 
energy levels next to loci, lowest is 1. 
1430 MHz 


= ee 8500 MHz R, and Ry operating points 
———— 10.7 KMHz . Be 
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with w ,/27=8500 MHz for the 1-3 transition and Ws/2% = 1430 MHz for 
the 2-3 transition with H=2100 6 and at 21° to the c-axis. 

With a ruby of a nominal ratio of Cr : Al=1 : 200 it was not possible 
to obtain an emissive condition, and all resonances saturated homogene- 
ously when applying the pump power. For Cr : Al=1 : 1000 emission was 
obtained at operating point R, whilst this was not the case at point R, 
for the same crystal. The experiments described in the following, were 
carried out on this last ruby. 


[3.9.1] Pump frequency 8500 MHz. Operating point Jie 
In fig. 3.26 the rise times of the signal, t,,, are plotted as a function of 
temperature. It is seen that they are close to 22/7’ ms. The temperature 
dependence appears to become less at the lower temperature range. 
For this operating point R, emission was observed. 


[3.9.2] Pump frequency 10.7 KMHz, operating point R&,. 
For the 1-2 transition at R, the absorption observed was much stronger 
than for the 2-3 transition at Ry, in agreement with a calculation using 
the matrix elements. 

However, on applying the saturating field, only a brief pulse of emission, 
lasting a few milliseconds was obtained. This can be seen in figure 3.20b 
where the 10700 MHz oscillator is pulsed on for 13 ms. every 40 ms., as 
indicated on the lower trace. The upper trace shows the signal reflected 
from the cavity at 1430 MHz, emission from the (unpolished) ruby causing 
a decrease in the reflection (downwards deflection) for the undercoupled 
cavity used. When the 10.7 KMHz oscillator is switched on, there is a 
short pulse of emission, followed by an increase in the absorption to above 
its equilibrium value. This was observed at both 1.4° K and 20° K. The 
time constant has been plotted in fig. 3.26. Other measurements then 
showed that subsequently the absorption decreased slowly, and after 
about 1 s. (at 1.4° K) it had fallen to a steady value close to half the 
equilibrium absorption. Even with high saturating powers (> 200 mW) 
the steady state absorption could not be further decreased. 

As the population difference only decreased to half the equilibrium 
value at 1430 MHz one expects absorption at the 2-3 transition (9270 
MHz), the idling frequency, to become negative when the pump is fully 
saturating. We therefore carried out a triple frequency experiment and 
observed the signals at 1430 MHz and 9270 MHz simultaneously. The 
2-3 signal was seen to decrease to a very small value, but not to turn 
negative, indicating that the pump transition was not fully saturated 
which is surprising because of the large pump power used (> 200 mW). 
This spike behaviour may be related to the intermittent oscillations 
observed for continuously pumped ruby maser at 9000 MHz by Maxkuov, 
e.a. [37]. 

The occurrence of the spike is confined to a region of the field and 
its angle narrower than in which the transferred modulation is visible, 
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° 
irae 5 10 20 30K 
Fig. 3.26. Characteristic times of signal absorption for ruby 


A Tor Pump 8.5 KMHz 
O Tm pump 10.7 KMHz 
O decay time spike (pump 10.7 KMHz) 


Also the spike decay time is a sensitive function of field and angle, being 
a minimum exactly at the operating point and increasing by a factor 
four, 2.5° further from the axis in which direction it became equal to T_ 
and where the spike started to disappear. Rotating H towards the axis 
made the spike disappear without changing the decay time. Increasing 
H from the value at R, on, increased the spike decay time; decreasing H 
from R, made the spike disappear. . 

The characteristic time t,,¢ of the signal at 1430 MHz is plotted in 
fig. 3.26 as a function of temperature in the range from 20.4° K to 1.4° K. 
One observes t,, in that range to be a constant above 4° K, while below 
4°K a T~- law is followed. 


[8.10] Discussion of the pulse measurements 

[3.10.1] Discussion of the results on KsCr(CN)g/KsCo(CN)g. 
According to the considerations in [1.6.2] one expects to find three time 
constants for sample 4, Only one was found as is evident from fig. 3.21. 
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The deviations from the straight line in that figure are due to a time 
constant of about 80 ms. in the apparatus. We will now relate the 
relaxation rates to the rise times observed by assuming that the fourth 
level does not take part in the relaxation processes. The time dependent 
equations (1.64) and (1.62) are then easily solved (ef. (1.65)) and one finds 
for the reciprocal of the two time constants left (1-2 signal and 1-3 pump 
transition) 


(3.05) A,2= Vie+Vis+Vos+ V Vie? + Via2+ Vos2+ Vie Vis + Vie Vos + Vig Vos. 
This expression can be simplified when it is allowed to expand the root 
and one then finds with r=Vj2/Vo3 and R= V43/Vog 


_ pg. Rr+rtkh 
\ a= Vo3 50 Rt4) 


(3.06) 
2Vo3(1+ R+r) Jie 


I 


lt 
For the time behaviour after shutting off the pump power we can try 
to identify the reciprocal rise time observed with 4; or with Ag. Inserting 
the values found by steady state measurements at 1.4° K for Uy, U3 
in [8.8.3] and 4=103/44 s-! for operating point I of sample A, results 
in a negative value of R when 4), and R=3.1 when Az is the reciprocal rise 
time observed. 
Only the last solution is acceptable, so we find: 


U12=0.6 8-1; Uo3=2.6 8-1; Uig=8.1 8-1 


and J,=1.2 s-1 for 7 =1.4° K. The values of U when plotted as a function 
of w on a double logarithmic scale are on a straight line so that U « w?-8. 

This proportionality should not be taken too seriously as according 
to (1.51) the matrix elements of the spin operators for the different 
transitions do enter. The expected frequency dependence for equal matrix 
elements is w? according to (1.53) and (1.60). From these experiments 
no information is obtained about relaxation rates involving the fourth level. 

During the pump pulse, Jz is large, while 2; approaches the limit 


= which equals 1.6s—1. The constant 1//1,,,,. characterizes 


the time it takes to establish a stationary population of level 2, the pump 
transition being saturated. This time constant was observed to be really 
so long, as it took roughly one second for the maser to start amplifying 
after stepwise applying the pump power. 

The reason for working with the particular concentrations used was, 
that for sample A maser action was obtained while this was not the case 
for sample B (cf. table X). This indicates that for A the spin lattice 
relaxation rates determine the individual level populations while in 
sample B the spin system comes through cross relaxations into equilibrium 
in a time shorter than a critical value determined by the spin lattice 


relaxation times. 
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When cross relaxations become important two situations may arise. 
In the first U> Uti, (the largest Uy present at that temperature), 
the time behaviour of the population difference of the signal transition 
is essentially determined by Uyj,,,. The second situation may be met 
with when the cross relaxation rate forms a bottleneck for the energy 
transfer, but is still fast enough to prevent maser action. So when we have 
(see [38]) 


(3.07) U.. on > Um om — Un ox 
where @x1=@s, and 
(3.08) Tig < Use) on 


As the Ujy’s are temperature dependent, while this is not the case for 
U,,, these conditions may be fulfilled at one temperature, but at a lower 
temperature the first situation, or at a higher one maser action may 
be expected. We suggest that in the measured temperature range below 
4° K the second situation exists for sample B with U, & 108/31 s7. 

Above this temperature an inverted difference in populations of transi- 
tion 1-2 or at least a non-equilibrium situation in the spin system seems 
also possible for sample B. At temperatures much lower than the ones 
measured at tT is expected to increase again (situation I). Condition 
(3.08) can be fulfilled by assuming that at least one of the relaxation rates 
from the fourth level towards the three levels just considered is fast. Theo- 
retically this is expected to be Uy4. 


[3.10.2] Discussion of the results on the ruby. The charac- 
teristic times Tt, observed at the operating point R,; (pump 10.7 KMHz) 
can be explained in the following way. Relation (3.07) holds for the whole 
temperature range measured, but below 4° K we have U,.>Uy,,,. AS 
Vijnax & Z’ condition (3.08) may be fulfilled at a higher temperature and 
this appears to be the case above 4° K, with U, » 1038/5 s-. 

As for operating point Rj, emission was observed, the inequality (3.07) 
does not hold. Cross relaxations appear to become important in the lower 
temperature range. The spike behaviour is not quite understood. If only 
three levels are concerned, relaxation processes will not increase the 
absorption above its equilibrium value. However, this could be explained if 
there were cross relaxations involving the fourth level. Since Teptce SL 8 
relaxation rate towards this fourth level would then have to be the 
determining factor. 


(To be continued) 
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PHYSICS 


ON THE PROBLEM OF NEUTRON DIFFUSION 
BY 


N. G. VAN KAMPEN 


(Communicated by Prof. H. FrEuDENTHAL at the meeting of November 28, 1959) 


1. Introduction 

The transfer of radiation through an atmosphere and the diffusion of 
particles (e.g. neutrons) through a medium are governed by the same 
transport equation. For convenience we shall use the language of particles. 
If the medium is a flat layer of homogeneous material, and if the individual 
scattering processes are isotropic and elastic, the general equation reduces 
in the stationary case to an integro-differential equation in two variables, 
viZ., 


+1 
(1) pA) + fe, u) = 4 | few’) du’. 


Here rt is the coordinate perpendicular to the surface of the layer, measured 
in appropriate units (“optical depth’); we put t=0 on the left surface 
and t=a on the right one. « is the cosine of the angle between the velocity 
of the particle and the t-axis and varies between +1 and —1; f(r, ) 
du is the density at the depth t of those particles that move in a direction 
between yw and w+du. Equation (1) is a balance equation; the first term 
gives the variation of the number of particles due to convection, the 
second term describes the loss due to scattering into other directions, 
and the third term represents the gain due to scattering from other 
directions into the direction p. 


Y 
es L 
2 
fo Gur jG 
T=0 Tzq 
‘a 


Fig. 1. The flat medium. 


When the factor j in the first term of (1) is replaced with 1, the resulting 
equation is of a very simple type, known from the theory of Markoff 
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processes '). It can easily be solved, subject to the boundary condition 
that the angular distribution /(0, w) of the particles incident on the left 
surface of the medium is prescribed. The presence of the factor Wes 
however, materially changes the structure of the equation, owing to the 
fact that ~ takes both positive and negative values. The correctly set 
boundary conditions are now quite different, as is intuitively clear from 
the physical interpretation of (1). It is now natural to prescribe the 
angular distribution of the incident particles on both surfaces of the 
medium, i.e., 


(2a) (0, u)=gi(u) for O<p<1, 
(2b) f(a, 4) =ge(u) for —1<p<0, 


where g; and gz are given non-negative functions 2). 

A simple, but important, special case obtains for very thick layers, 
a —>oo. It turns out that in this limit the detailed distribution of the 
incident particles on the surface t=a is immaterial; only one parameter 
remains in the solution, namely the total current. This case of a semi- 
infinite medium is the only one for which exact and explicit solutions 
have been found. If, moreover, there are no incident particles on the left 
surface, gi(u)=0, but only a given total current, one has the simplest 
possible case, known as Milne’s problem. It has been reduced by Milne 
to a one-variable integral equation, which was solved by Wiener and 
Hopf. An explicit solution for the semi-infinite medium with arbitrary 
incident particles has been constructed by Chandrasekhar. 

These solutions, however, are obtained by methods specially designed 
for the present problem and have no apparent connection with the more 
familiar methods of mathematical physics. The object of the present 
paper is to furnish a more straightforward treatment of equation (1). 
First a general expression for the solution of (1) is obtained in the form 
of definite integrals involving two arbitrary functions (section 3). These 
arbitrary functions have to be determined from the boundary conditions. 
This is carried out for the case of a semi-infinite medium in sections 4 to 6. 
The one artifice used is the ‘‘Ansatz’’ (13). The mathematics used in 
adjusting the general solution to the specific boundary conditions is akin 
to the Wiener—Hopf technique, or rather to its more general formulation 


1) H. Onsen [Kong. Norske Vid. Selsk. Forhandlinger 28, no. 2 (1955)] has 
pointed out that this approximation is the basis of the various small angle treatments 
of multiple scattering of fast particles; see S. Goupsmir and J. L. SAUNDERSON, 
Phys. Rev. 57, 24 (1940); G. Moimrp, Zeits. Naturf. 10a, 177 (1955); H. 8. SnyDER 
and W. T. Scort, Phys. Rev. 76, 220 (1949). 

2) It is not the purpose of this paper to study the mathematical problern of 
proving the existence and uniqueness of a solution of (1) with boundary conditions 
(2), Accordingly we do not try to specify the exact restrictions to be imposed on 
gi en ge, nor those to be imposed on the functions 41 and y2, which appear later 
in the text. 
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by Hilbert, but it will not be necessary to employ the general theory. 
Although the results obtained are not new, it is felt that a more streamlined 
derivation of old results may be useful for future generalizations *). 


§2. Preliminary remarks 
The net current of particles through a plane at the depth t is 


41 
(3) J= § uf(t, udp. 


By integrating (1) with respect to mw one finds that J is a constant 
(independent of r), as is physically obvious. The transport of momentum is 


+1 


(4) K(t) = J pf (t, w) dy. 


One easily finds from (1) that dK/dr=—J, so that 
5) K(r) = — J(t+1T1). 


The integration constant t; depends on the boundary conditions (2) and 
can only be determined when the complete solution has been found. 
In particular, in the Milne case 1 can only depend on J; but as K is 
necessarily linear in J, t; must be a universal constant, which we denote 
by t. It determines the second moment of uw for the emerging particle 
distribution in the Milne case. 

The fact that + does not occur explicitly in (1) suggests that one try 


(6) f(t, uw) =e" u(y). 
Substitution in (1) yields 
es 
ez , , _ const. 
(7) u(u) = y=, J wu’) dw ia 


This will constitute a solution, provided that, firstly, the parameter s 
is such that the denominator does not vanish for —l<p<+1; and 
secondly, that the following consistency condition is satisfied, which is 
obtained by integrating (7) with respect to y, 


+1 

a dues l+s 

(8) 1=4f Ao = glo Se. 
Si 


Each root of this equation gives rise to a solution of the type (6), (7). 
Unfortunately there are only two roots, coinciding at s—0. They give 


‘) In fact, the method can be extended to certain types of anisotropic scattering, 
including Rayleigh scattering. 
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rise to two somewhat trivial solutions, viz., 


(9) He.) = [7 ]_ = 
and 
(10) fee) = [I] = eee 


The most general solution found in this way is a linear combination of 
both with two arbitrary constants 1), 


(11) f(t, w) = A—Biu—2). 
One readily checks 
(12) J=—-—$B, K=}#A+Br), 1=A/B. 


Obviously the two constants are not by far enough to satisfy the boundary 
conditions (2). In the next section the Ansatz (6) will be modified so as 
to obtain the general solution of (1). It will appear, however, that (9) 
and (10) are the only solutions that contribute to J and K. 


§ 3. The general solution of eq. (1) 
As a generalization of (6), (7) we make the Ansatz 


1 e—st 


201 J 1l—syu 
W 


(13) f(t, u) = p(s) ds. 
where the function ¢(s) and the integration path W are to be determined 2). 
Substitution of (13) in (1) ee 


(14) J e-"[1= =, log — = | v8) 8) ds = 0. 
Ww 


This condition can be fulfilled by taking for W a closed contour and 
choosing ¢(s) such that [ ]y(s) is a holomorphic function in the interior 
of W. Of course, ¢(s) itself must not be holomorphic inside W, because 
that would give f(z, ~)=0. Each zero of [ ] yields a possibility of meeting 
these requirements, by choosing for g(s) a function with a pole coinciding 
with the zero of [ ]. We use the abbreviation 


(15) = log 572-1 = Als). 


l—e 
As mentioned before, A(s) has a double zero at s=0: 
(16) A(s) = 482+ O(s4) 

Hence a solution is obtained by putting 

(17) p(s) = As = Bs, 


1) Since f cannot be negative, one must have 4 > > |B|+4a(|B|— 
2) It should be noted, however, that W must not contain any ae on the 


real axis with s >1 or s<—1. 
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and taking for W a circle about the origin. This leads to the solution 
(11) of the previous section. In order to find more solutions it will be 
necessary to study the analytic behavior of A(s). 

A(s) is holomorphic in the complex s-plane with the following two 
cuts along the real axis, 


(18) er pat Be Ce:8s<—l. 


On the upper edge of Ci, A takes the values 


(19) A(s+) = 5; (log + mi) — 1. (s>1) 
iho 
Oo " 
CW SLLLLITIL, GIS LO ALLS 
2 Garey ras Zo LSS =1 MLL LLL Cy 
Wo og—ie€ | 


Fig. 2. The complex s-plane. 


The values on the lower edge follow from A(s*)=A(s)*; consequently 
there is a jump on C; given by 


(20) A(s+)—A(s—) = + at/s. (s>1) 


The behaviour on C2 can be found from this by using A(—s)= A(s). The 
behaviour at infinity is sufficiently described by 
(21) lim A(s) = —1. 
|s|—+0o 

We shall make the argument of A(s) unique by taking it zero for 
0<s<1. One then finds that it increases from 0 to a along the upper 
edge of Ci, going from s=1 to s= +oo, Hence the argument of A(s) 
tends to +a in the upper half plane, and to —z in the lower half plane. 
This remark will be of future use, and also shows that A(s) has no other 
zeros than the double zero at s=01). 

Now take for p(s) a function that is holomorphic in some neighbourhood 
O; of Ci, but has a jump on (C satisfying 


(22) p(s + )/p(s—) = A(s—)/A(s +). 


Clearly A(s)p(s) is holomorphic in O;, including C;. Hence, if one takes 
for W the loop W, enclosing 0, the condition (14) is fulfilled without 
(13) vanishing identically. An obvious solution of (22) is 


(23) p(s) = xa(s)/A(s) 


1) Consider the closed contour formed by W; and We and two large semicircles 


(lig. 2). The argument of A(s) increases by 47 along this contour. Hence there are 
just two zeros inside. 
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where y1(s) is holomorphic in Qj, including ©). Thus we have found the 
following class of solutions of (1), 


24 : eee ae) 
( ) I(r, ph) = Oni | l—su A(s) ds. 


1 


Here yi(s) is an arbitrary function, holomorphic in some neighbourhood 
of Gye 

Actually, this restriction on y;(s) is immaterial, as the neighbourhood 
may be arbitrarily narrow. To make this evident one may collapse Wy 
onto C; so that (24) becomes for —1<y<0, 


) — 8T x1(8 
( >} -4/z l—syu sR(s 
ve 
where we used the abbreviation 
* La s+l1 m2 
(26) R(s) = A(s+) A(s—) = (5 log St 1) +i. 


For 0<y<1 there is a pole on C; at s=1/u, which gives rise to an 
additional term: 


— AQ/u,+)+A(1/e, — — om HO 
(27) f(z, m) = ep na(z)¢ oe SS sR(s 


In the integral the principal value has to be taken at s=1/y. 
Similarly the cut Cz furnishes a second class of solutions 
1 e—s (8 
(28) f(t, m) = 55 ii ear ee 
Ws 
with another arbitrary function y2(s). The most general solution that has 
now been found is a superposition of (11), (24) and (28), with arbitrary 
A, B, m1, y2. A direct calculation shows that both J and K are zero for 
the solutions (24) and (28); hence the transfer of particles and momentum 
through the medium are still given by (12). If the medium is very thick 
compared with the mean free path (a > 1), the connection between the 
arbitrary elements A, B, 71, 72 in the general solution and the boundary 
conditions can be understood in the following way. (24) vanishes at least 
like e~* and is therefore only required to describe the situation near the 
left surface. Similarly (28) is confined to the vicinity of the right surface. 
Between those two boundary layers (11) is a good approximation. Each 
boundary condition determines one function x and in addition gives rise 
to one linear relation between A and B. Hence both boundary regions 
are necessary to find the constants A and B, which determine the solution 
in the bulk of the medium and the transfer of particles and of momentum . 
To prove that the solution obtained above is the most general one, 
it has to be shown that A, B, 71, v2 can be adjusted to give any set of 
prescribed distributions (2) for the incident particles. This will be done 


7 Series B 
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in the remaining sections, but only for the case of a semi-infinite medium, 
q—co. Then the term (28), which refers to the right surface, is immaterial 
at any finite depth t, and must be omitted. The connection between 
A and B arising from the boundary condition (2b) must now be replaced 
by an explicit relation. A suitable way to achieve this is to regard the 
current J as a given quantity. In the next section it is shown that with 
given g; and B and with y2=0, there is indeed just one solution for 71 
and A. 


$4. Adjusting the solution to the boundary data for a semi-infinite medium 


The general solution for a semi-infinite medium is 


est 


(29) H(t, uw) = 55 o(s)ds+ A+ B(r—y). 


2x1 J 1—sp 
Wy 


The boundary condition (2a) is 


1 ds 
(30) oni | a +A—Bu = gil), (O<p<1) 
Wy 


or, putting w=1/o and o-1!9\(0-!)=G(oa), 


1 ols)\ds A B_ 
Bh oer aver os Shae (o> 1) 
Wi 


(31) 


Clearly this equation requires G(c) to be holomorphic in Q;, but this is 
again no actual restriction on G, because QO, is arbitrarily narrow. Our 
task is to solve this equation with given G(a) and B. 

We shall employ the fact that the first term in (31) has the form of a 
Cauchy integral. For o inside the loop Wy one has the identity 


1 poles 1  ole+)—9le—) 
9(8) = — 55 | SPS + [SE as 
(32) pis s fey Tt : 8 oO 
= gi(o) + ui(e). 


Here g(a) appears as the sum of two parts. The first part, x(a), is 
holomorphic in Q;, including C;, but may have any number of singularities 
outside O,; the second part gu(c) is holomorphic everywhere but for a 
discontinuity on 0;. The first part is determined by (31), 


(33) gio) + 4 — 2 = Go). 
The second part must be so chosen that (22) holds, i.e., 
(34) A(s+) pu(s+)—A(s—)gn(s—) = —(ai/s)gx(s). (> 1) 


Thus the problem reduces to finding a function gp(s) whose discontinuity 
on C; satisfies the inhomogeneous equation (34) 


99 


This is the problem of Riemann—Hilbert, which can be solved by 
methods developed in the theory of singular integral equations. One first 
solves the homogeneous equation 


(35) A(s+) G(s +)— A(s—) B(s—) = 0, 
or 
(36) log P(s +) — log D(s—) = Log {A(s—)/A(s+)}, 


by putting (s not on (Cj) 


(37) log ®(s) = 55 | Log (Se \ =. 


This integral only exists if Log denotes that branch of the logarithm 
that vanishes at infinity. According to the above stated properties of 
A(s) this entails 


(38) lim Log {A(t—)/A(¢+)} = 2a. 


It follows that log ®(s) tends to infinity for s > 1, in the same way as 1) 
—log (s—1). Thus @(s) has a pole at s=1, which will presently become 
important. Otherwise @(s) is finite, different from zero, and tends to 
1 at infinity. 

With the use of this function ®(s), the inhomogeneous equation (34) 
can be solved for gyn(s). More directly, however, one may use Ps) to 
write the equation (22) for the total function g(s) in the form 


ger) o{s—) 


.) Diet) De)’ 


This shows that (s)/®(s) is continuous on Ci, so that one has, for s 
inside W,, 


y(s) g(c) 
(40) a= zn | Bidy 


Writing in the integral y(c) Dae one readily sees that the part 
with gr(o) vanishes, because gyr(c)/®(c) is holomorphic outside W,. 
Hence (40) determines g(s) in terms of 91: 


i gi(a) do ae 4 A B — do 
Me) Qa a) D(c) Ont | (4c) a 7 a) P(c)(o—s) 
Wy 


This completes the solution of the problem, but for one complication. 
As @(s) has a pole at s=1, it follows from (14) that the same is true for 
g(s). This would entail a pole of A(s)p(s), which is not permissible, because 
A(s)p(s) has to be holomorphic in order that (14) shall vanish. This pole 
can only be cancelled by causing the integral on the right of (41) to 


1) A more precise characterization of this singularity is given in the appendix. 
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vanish for s—1. That can be achieved by an appropriate choice of 4; 
thus the constant A is determined by the condition 


do 
W, 
To summarize the solution: First B is determined by the given total 
current J according to (12); next A is to be found from (42) where @ is 
given by (37); then ¢(s) follows from (41) and finally f(t, «) results from 
(29). In the next section various expressions for G are treated separately 
and for each of them the above equations are cast into a form that is 


more accessible to actual computations. 


§ 5. Milne’s problem 
No incident particles: @=0. The current J = —%B is given, A has to 
be determined from 


do tr da 
(43) A | oD(a)(c—1l) B | a? P(a)(¢—1)" 
Wi Wy 


By shifting the contour W, to the left, one finds for these integrals 
—@(0)-1 and —@(0)-!+ @'(0) @(0)-? respectively, so that 


_A_,_ (0) 


With the aid of (37) it is shown in the appendix that this can be written 
in the form 


ee) fee ene ae 
(45) ba aaa i) [; (1 +¢)(t —are ag - 
0 


This expression was given by Placzek and Seidel and computed to be 
0.71044609. In the same way one finds from (41), using (44), 


(46) Beg Nea 


Hence the solution of Milne’s case is 


B mee Jem, 
(47) f(t, 4) = ni B(0) orn ~=* P(s)ds+ A+ B(r—). 


Wy 


For t=0, w>0 the integration path may again be shifted to the left, 
with the result 


(0.4) = ge {®'(0)— (1) (0)| + A — Bus 


owing to (44) this vanishes, as it should. The emerging angular distribution 
is obtained from (47) for t=0, u<0; one finds, using (A11) 


(48) HO, #) = gop AH) (2) = —4VB IU + [uly 0(— a): 
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This shows that the function @ is essentially the angular distribution of 
emerging particles in Milne’s case (“law of darkening”). On putting 
J =—1, eq. (48) becomes identical with eq. (46) of Placzek and Seidel. 
For t>0 (inside the medium) the expression (47) must remain as an 
integral, but it can be made a real integral by collapsing W, onto Cy}. 
The general result was already given in (25) and (27); one merely has to 
substitute 

—31V37 


mls) = =~ —— A(s) Os) = (as) =e) 


were (A10) has been used, and to add the terms 


A+ B(t—p) = — 3d (t+ 7t0—p). 


§ 6. Arbitrary incident particle distribution 


Every solution of the problem for a semi-infinite medium can be regarded 
as a superposition of the solution of Milne’s problem with the prescribed 
value of J, and a solution with J=0 and the prescribed incident 
distribution g;. In this section we discuss the second solution; accordingly, 
we put henceforth B=0. (There is no difficulty in treating the general 
case with BAO in the same way, but the formulae would become 
unnecessarily lengthy.) 

If nothing is known concerning the analytic nature of gi(u) (e.g., if 
gi(u) is given numerically), the integral (41) cannot, of course, be 
simplified. The only thing to do is to reduce (41) and (42) to real integrals 
along C; (in the way of (25), (27)) and to evaluate them. On the other 
hand, if g;(~) is given in the form of an analytic expression some further 
reductions are possible. We shall briefly discuss two cases, viz., the case 
that gi(u) is a polynomial, and the case that gi(u) is a delta function. 
These two cases also provide approximation methods, respectively for 
smooth and for sharply peaked gi(y). 

First suppose that gi(u) is a polynomial 


(49) gi(t) = dcme™, G(c) = dS emo-™, 


The first integral in (41) becomes 


1 G(a) do 
(50) sai | Bm goa = Dom Pale) 
W, 
where 
if do 1 1 
(51) Yml8) = 57; J amti@(a)(o—s) amt Eas 
W: 


Here [1/®(s)]m indicates the polynomial of degree m consisting of the 
first m+1 terms of the power series expansion for 1/®(s). Hence the 
Wm(s) are polynomials of degree m+1 in s7}, whose coefficients involve 
the derivatives of ® at the origin. These can be calculated from the 
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definition of ® (see appendix). It turns out that each derivative of odd 
order requires the evaluation of a definite integral (A18), while the 
derivatives of even order can be found by means of the symmetry relation 
(A10). Actual calculations will therefore only be possible if the degree 
of the polynominal is not too high; but, in principle, this case can be 
treated as completely as the Milne case. 

Secondly, let gi(u) be a delta function, 


(52) gil“) = H(u—fo), (a) = o00(a — 90). 
In order to treat this case we have to replace 6(¢—o9) by an approximate 
expression which is holomorphic near C;. The most convenient one is 
00 é _ 00 1 i 1 
(63) Go) ~ 7 (a—oo)?+e2 sass ae ——;): 
where ¢ is a small positive number, which will eventually tend to zero. 
Then 
G(o)do 1 l 
A RD EH Eestighen BY ena ee Me ear Ae EW OS, Me le es 
(54) D(a) (a —s) 7 D(do +7e) (oo +7e — 8) D( oo —ie) (60 —1e —8) . 
Wy 
This can be used to find first the constant A, subsequently the general 


expression for f(t, “), and finally the emerging distribution /(0, u) for 
<0. The calculation is performed in the appendix with the result 


(55) f(0, w) = $u0(1+ Mo) o(- -) o( - va) ai 


Clearly this result makes it possible to write down the emerging distribu- 
tion for arbitrary gi(u) in the form 


1 
1+ 1 1 1 
(56) (0, #) =F @(— To) f BOE) w(— 2) gr (ua) do 
0 


However, owing to the difficulty of performing the integration, this 
formula will primarily be useful in the case of a sharply peaked incident 
distribution. 


APPENDIX 
The function D(s) 


Concerning the general features of ®(s), it follows from (37) that ®(s) 
is holomorphic and different from zero everywhere outside (1, and tends 
to 1 at infinity. For real s less than +1 it is real. Moreover, as the function 
—t Log {A(t—)/A(t+)} vanishes for t= +00, and has a negative derivative, 
it is always positive; hence @(s)>1 for —co<s<l. Moreover, the 
imaginary part of log @(s) is positive in the upper half plane and negative 
in the lower half plane, so that ®(s) increases monotonely from +1 
to +oo as s moves along the real axis from —oo to +1. 
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S=0 S-| 


Fig. 3. The function ®(s) sketched for real s. 


Next we write (37) as a contour integral. Let log A(s) denote that 
branch of the logarithm that is holomorphic in the right half of the 
s-plane excepting Ci, and real on 0<s<1. Then for s>1 


(A 1) log A(s+) —log A(s—) = Log 7 a+ 271. 


As log (1—s) has a jump [log (1—s)]* = —2z7 on (Ci, 


(A 2) [log {(1—2) A(s)}]+ = — Log =. 


We insert in { } a factor (1+3s)/s?, so that the logarithm is holomorphic 
at s=0 and tends to zero at infinity. Thus (37) may be written in the form 


ao—8s 


(A 3) log ® ~ a | 18 [A (c) Baas 


Here s lies outside W,; on moving s across W, to the inside one gets a 
separate term from the residue: 


(A4) log ®(a) = — log (25* A(6)|— 4, f log [~S* A(o)| =. 


The integral is holomorphic in 0;, so that the singularity of @(s) at s=1 
is specified by 


(A 5) (1—s) A(s) &(s) is holomorphic at s=1. 


This justifies the somewhat loose argument at the end of § 4. Indeed, 
after A has been solved from (42), the right-hand side of (41) contains 
a factor (1—s). Hence (41) states 


(A 6) p(s) = ®(s) (1—s) x holomorphic function. 


Because of (A 5) this implies that A(s)g(s) is holomorphic in O;, which 
is the necessary requirement in order that (14) shall vanish, 
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Equations (A 3) and (A 4) may be looked at in the following way. 
Let the function 
1—s? 
82 


(A 7) L(s) = — log { A(s)}, 


be decomposed, similarly to (32), into a term Jy that is holomorphic in 
O;, anda term Ly that is holomorphic everywhere but for a discontinuity 
on C;. Then (A 3) and (A 4) state that 


(A 8) log ®(s) = Im(s) = L(s) —Lx(s). 


As L(s) is even and has no other singularities than the jumps on C; and 
Cs, one has 


(A 9) In(—s) = In(s). 
Hence 
(A 10) (8) D(—s) = eh = ~ 


In particular, by taking s=0, one finds from this 
(A 11) (0) = y3, 
the sign being determined by the fact that ®(s)>1 for s>1. 
The relation (A 10) also yields the following important property of ®(s) 


9 1 a 1 . tse 
ogee! Dist) Oe) “ s 


P(—s), (s>1) 
which will be used in (A 20), (A 22), and (A 25). 


The polynomials wm(s) 


In order to calculate the polynomials ym(s) defined in (51) we write 
the expansion of 1/®(s) =e) 


(A 13) ee 7 (1 —sLan' + 5 n'*— Ln"} — F [Lan — Zar’ Lar” + Ln") 


4 
+ $ [a's Lay"? Dan" + 4Don' Inn" + 8L1"2— Ln”) +...) 


The derivatives of Iq of even order can be found by expanding both 
members of the identity 


(A 14) In(s) + Im(—s)= —log (1—s?)—log {s-2A(s)}, 
with the result 


(A 15 Deh iOl wee we )_ 22 
) n’/ B Dy""""/4! Tee 
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To find the derivatives of odd order, the explicit expression (A 3) for 
® must be employed. One has for m=1, 3, 5, ... 


Ly”) = [(=)” log (a) ]_ = - | tog (= A(o)| 


8=0 270 


_ _. (m1)! —26 2. A'(a))\do 
(A 16) = | eet ci 
Wy, 


(m = 1)! — idan = i (1 —o?) = 


A(o) gm-1° 


I 


The best way to make the remaining complex integral accessible to 
numerical calculations is by shifting W , onto the imaginary axis. The 
multiple pole at the origin can be made harmless by subtracting a number 
of terms, which do not contribute to the integral. Thus for the integral 
n (A 16) we write 


1 o = do 
(A 17) sa) goa Oe a | gmt? 


where the coefficients « are so chosen that [ ] has no terms of lower 
order than o”+1,. Then this integral equals 


nr t—larc tant — 1) 


—1)#(m-1) t2 
(A 18) ee | (Ge Fo oe” ae | aa 
0 


Once these integrals are evaluated numerically, the successive derivatives 
@(™)(0) can be found. In particular, for m=1 the result is given in (45). 

Denoting the integrals (A 18) by 1—Jm one finds for the expansion 
of 1/®(s) =e) 


3-4(1—lis+(3h2—z)s*— (Gh? — 5 h+gla)s* + 
(A 19) i ; 
+ (Ght-pht+5 =I, [3 — aaa) 8 4 ies, 


from which the successive pm(s) follow by means of (51). 


The case determined by (52) 
By taking s=1 in (54) one finds with the aid of (42), (A 11) and (A 12) 


Qi 1 1 Lees oo+1 = 
(A 20) — 74 ~ ay = a Baas | i we O02 P(— 00). 


Substituting (54) and (A 18) in the general expression for f(t, “) 
il 00 
(A 21) f(t,p) = 4 + 5 i (6) a0 = 1] se 


Wi 
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where {} is the same expression that occurs in (54). For t=0 the 
integration path may be shifted to the left, yielding for “<9 


10.) = =°(;) + 33,59 (7) (aaa eteSIn ~ 


Bae eR TI 
: uh Peed erie) 


- tp — 00) &(=) = sae 


which is the result (55). It can readily be checked that this gives indeed 
an emerging current of the same strength as the incident current corre- 
sponding to (52): 


(A 23) Tout = J u f(0,u) du = = po. 


(A 22) < 


In this derivation it was assumed that o941, and, in fact, this value 
had to be excluded already in (52), as d(u—1) has no clear meaning. 
However, it may now be noted a posteriori that the final formula (55) 
is perfectly regular for wo=1l. Thus, if d(u—1) is interpreted as 
d(u—1+e’), the formula (55) still holds. Consequently, the angular 
distribution of the reflected particles, in the case of an incident beam of 
unit strength perpendicular to the surface, is 


(A 24) f(0,u) = ®(—1) Bu). (u“< 0) 


Comparison with other work 


Chandrasekhar’s function H(u) is related to our @(s) by H(u)= 
=(1+y) ®&(—1/u). His integral equation for H(u) can be derived by 
applying Cauchy’s theorem to the function s{(1—s) ®(s)]-!. This function 
tends to —1 at infinity and has a discontinuity determined by (A 12); 
hence 


(A 25) mae {Ss or ol. 
0 


Substituting s= —1/u and s’=1/y' one obtains 


A(u') a 
Ba 1 | a 


Another relation is obtained by applying Cauchy’s integral theorem to 
[(1—s) ®(s)]-, viz., 
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also given by Chandrasekhar. The same method applied to 
[sm+(1—s) D(s)] 


yields a relation connecting the polynomials y» with the H-function 


1 
( arr eae) 
fi 4 PSH ae  $ va(—2) 

0 

An infinity of similar equations may be obtained in this van from (A re 
Placzek and Seidel used the functions t4(s), r-(s) and 1(s)=1+(s)/t_(s 

They correspond to our e-Z1), eZ), and e-£) respectively, so ne 
t_(s) =@(—s), and t+(s)=[®(s)]}-. 
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PHYSICS 


ON THE CAPACITY OF A SPHERE INSIDE A CUBE 
BY 


F, W. DE WETTE 


(Communicated by Prof. L. van Hove at the meeting of November 28, 1959) 


Summary 


A condensor consisting of a conducting sphere placed concentrically inside a 
conducting cube is considered. A rapidly converging expression for its capacity C 
is derived with the aid of a method which has been developed elsewhere for the 
description of electrostatic fields in ionic crystals. Two approximations for C are 
discussed. 


The study of electrostatic fields in ionic crystals leads, as a by-product, 
to the solution of an apparently quite different problem, viz. that of the 
electrostatic field in a condensor consisting of a conducting sphere placed 
concentrically inside a conducting cube. A direct solution of the Laplace 
equation inside the field region, 7.e. the region between the sphere and 
the cube, is hampered by the geometric dissimilarity of the two boundaries. 
A solution can be found, however, when one recognizes the equivalence 
of this problem with that of determining the potential inside a cubic 
lattice consisting of conducting spheres, with equal radii and equal net 
electric charges, the signs of these charges being alternately positive and 
negative (see fig. 1). The field inside the condensor is identical with that 
inside an “elementary” cell of the cubic lattice (i.e. a cube containing 
one sphere, indicated with solid lines in the figure) 1), because in the 
condensor as well as in this elementary cell of the lattice, the surfaces 
of the sphere and the cube are equipotentials. 

The problem of determining the potential field inside the lattice can 
be solved with a self-consistent procedure. As a result of mutual inductions, 
the charge distributions on the spheres will not be spherically symmetric 
but they will have the cubic symmetry of the lattice. The potential 
field due to such a non-spherical charge distribution can be described as 
arising from a series of higher order multipoles residing in the center of 
the sphere. The potential fields due to all the charge distributions in the 
lattice, or, as we can now say, due to the net charge and the higher order 


') Since the elementary cell contains only one sphere, we have to distinguish 
between positive and negative elementary cells, containing a positively charged 
sphere and a negatively charged sphere respectively. This elementary cell is not 


a unit cell of the lattice of charged spheres. The latter contains at least one positive 
and one negative sphere. 
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multipoles at all the lattice sites of the lattice, add up to form the potential 
around any single sphere in the lattice. This potential must, at the surface 
of such a sphere, exactly induce the charge distribution which is present 
there. J.e. it must induce on the sphere just this same set of multipoles. 
With the aid of this condition a set of self-consistent equations can be 
derived which determine the higher order multipoles. Once these multipoles 
are known, the potential in the lattice is completely determined. 


Ve +Vo ) (V2 Vo, 
| | ee ee 
OO!@ 
| | : 


; i) 
woO@O} 


HoO@ Be 


\V=07% 


Fig. 1. Right: The condensor consisting of a sphere placed concentrically inside 

a cube. Left: A cubic lattice of charged conducting spheres. The straight lines 

represent equipotential planes of zero potential. The elementary cell is drawn with 
solid lines. 


This self-consistent method has been developed elsewhere [1, 2] for 
the more general problem of the determination of the electrostatic fields 
in arbitrary ionic lattices and for a subsequent evaluation of some internal- 
polarization effects in such crystals. For the present purpose we will only 
briefly indicate the derivation. 

The potential inside the condensor (or the elementary cell of the cubic 

lattice) may be expressed as 
(1) V(R) = & (Geert Ann B) YEW 9). 
R(R, 8, y) indicates the fieldpoint and Y;,, is the normalized spherical 
harmonic of degree k. The first term in the righthand side of (1) describes 
the potential due to the charge distribution on the surface of the sphere. 
Qr, n is a multipole moment of order 2. For a charge distribution on a 
sphere (radius r) with density o(9,¢) Qx,n is defined by 


(2) Quan = seat || (0, 4) r+ Yo. n(O, 4) sin 0.d0 dg. 


(r(r, 0, ¢) gives the position on the sphere). The second term in the 
right-hand side of (1) gives the potential due to the charge distribution 
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on the inside of the cube (which is identical with the potential of all the 
other charged spheres in the lattice). It can be shown that Ax, n has the 
following form (cf. [1, 3]) 


4a ] , n 
\ A... = Seer * grt SenlO[ky, 2) + 


1 1 , 
/ 9 > La a qkt+lF1 Qi, Sé41.n—m(0|k;, 3). 
1>0,m 
a is the length of the cube edge (i.e. the lattice distance in the corre- 
sponding cubic lattice). The quantities S’ are lattice sums for simple 
cubic lattices. They are defined as follows (cf. [1, 4]) 


1 
(4) Si,,r(O| ky, +)= {=I etnies oti Y 5 (9,, a). 


Ay, Ae, Ag = — 00 


2 stands for the three indices 2), 22 and A3 which can take all integral 
values between —co and +o. The prime on S and 2X means that the 
term A, =A2=43=0 is excluded from the summation. e,(0,, 9,, ¢,) is the 
lattice vector. Fi, m:; x, n is a numerical factor defined by (ef. [1, 3]) 


—(—])ems ; 21+1 _ (k+l+n—m)! (k+l—n+m)! 
(8) 2imina=t— 1) 2\/x (2k+1) (2k+21+1) (+m)! (I—m)! (k-+n)! (k—n)!" 


The lattice sums of the type (4) were first introduced in reference [4], 
where a rapidly converging method for their numerical evaluation has 
been developed. Values for various of these lattice sums for cubic lattices, 
calculated with this method, are listed in [1, 2]. It is the availability of 
these numbers which makes the present treatment of the condensor 
problem possible. 

The condition that the potential V be constant on the surface of the 
sphere, viz. V(r)=Vo, means firstly that the coefficients of the terms 
with k>0 in (1) must vanish. This leads directly to the self-consistent 
equations 


(6) Qn =— rt A,,. (4=1, 2, ...) 


Secondly, it yields an expression for Vo in terms of the net charge and 
the higher multipoles, viz. 


1 /Qo,0 1 l 
wf Se eee | et = = ——— 
(7) Vo Al Pe + Apo) , Te Ayo, 
where, for simplicity, we have chosen the net charge (422)*Qo,0 to be 
unity. Since the potential of the cube is zero when it is grounded 
(for the lattice this follows from the symmetry of the lattice arrangement, 
cf. fig. 1) we have for the capacity C of the condensor 


(8) GC=——__charge 
potential difference Vo" 


Walal 
From (3), (7) and (8) we find 


1 N= 
(9) C= )r4+ es eta Oe ta (Ok, 40 
where we used the fact that V/420 86, o(0|k;, 4) is just the well-known 
Madelung-constant xy (ef. [1, 4]). For a simple cubic lattice «4 —1.7476. 

The self-consistent equations (6) describe the induction of an infinite 
set of higher multipoles. The importance of these multipoles decreases 
rapidly with increasing k and this causes the infinite series in (9) to 
converge. In actual calculations the equations (6) are therefore cut off 
at some low value of k. Moreover, due to the cubic symmetry of the 
electric potential inside the condensor, only certain k, n-combinations 
give multipole moments different from zero. These combinations are 
k=4, n=0, +4; k=6, n=0, + 4; k=8, n=0, +4, +8; ete. The 
capacity can thus be given by various approximate expressions depending 
on the k-value at which the self-consistent equations have been cut off. 


1. The monopole approximation is obtained by cutting off after k=0. 
This assumes the charge distribution on the sphere to be spherically 
symmetric. This approximation is fairly good due to the fact that the 
first non-zero multipole moment occurs for k= 4. We find from (9) without 
further calculation 


_ 1, ow 
(10) Onon = + 


—1 1 1.7476 )-1 
r a 


2. The hexadecapole approximation is obtained by cutting off after 
k=4. Solution of the self-consistent equations and substitution into (9) 
gives 


; 1 1.7476 16.468 = 
(11) Cy a ( r , 


yy a(2)— 234.034 


The third term on the right of (11) is the hexadecapole contribution. 
Better approximations would add more terms to (11), but this would 
require knowledge of the values of lattice sums which are not readily 
available, although they could be computed if necessary without too 
much labour (cf. [1, 4]). We now roughly estimate the convergence of the 
successive approximations for ( from (10) and (11) only. It is obvious 
that the convergence will be worst when all multipoles have maximum 
values. This happens when the sphere is as large as possible, 7.e. for 
a/r=2. In that case we find from (11) Cigy*=5.1797 a, whereas from 


hex 


(10) 4C2 > = 3.96130; Hence (CR - CF Cree = 2” %, obviously the 


mon 
accuracy of (10) is very poor in this case. The percentage deviation of 


C-, from the real value of C, which is caused by the neglection of all 
multipoles with k>4, may very well be of this same order of magnitude. 
However, the convergence improves rapidly with decreasing r. For 


a/r=2.4 (11) gives C,,,= 1.5489 a, while from (10) Cyyo,= 1.5328 a. We 
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now have a difference of just 1 °%. For a/r = 3.0 these numbers are 0.7990 a 
and 0.7984 a, respectively, showing a difference of only 0.08 %. It seems 
reasonable to expect that the relative correction due to the higher 
multipoles (k>4) will never exceed the ratio (Cye,—Cmon)/Cmon- 

We conclude by pointing out the similarity of (10) and (11) with the 
expression for the capacity of a condensor consisting of two concentric 
spheres, viz. C =(1/r—1/R)— (where r and R are the radii of the inner 
and outer sphere respectively). In fact (10) may be written as 
Cnon = (1/7 —1/Reg)-+ where R,,=0.5722 a is the radius of an effective 
sphere replacing the cube. R,, lies between the radii of the inscribed and 
circumscribed spheres of the cube, as was to be expected. Also, the 
appearance of the Madelung constant in these formulae is not surprising 
if one realizes that application of the method of image charges to the 
condensor leads to just the lattice of alternately positive and negative 
conducting spheres, which we used in our description. 

The author wishes to thank prof. B. R. A. Niyboer for pointing out the 
relation of this condensor problem to the crystalline field problems which 
were treated elsewhere. 
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